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SUMMARY 


A  FORTRAN  IV  computer  program  to  generate  the  geometrical  input  for  the 
TACTl  program,  which  calculates  transient  and  steady- state  temperatures,  pres¬ 
sures,  and  flows  in  a  cooled  turbine  blade  or  vane  with  an  impingement  insert,  has 
been  developed  and  is  described  in  this  report.  The  spline  curve-fitting  technique 
used  in  previous  NASA  codes  to  calculate  the  flow  field  around  an  external  blade  shape 
is  used  in  the  present  program,  TACTGRID,  in  order  to  generate  the  internal  geome¬ 
try  of  an  impingement-cooled  blade. 

TACTGRID  constructs  the  blade  internal  geometry  from  the  previously  specified 
external  blade-surface  points  and  newly  selected  wall  and  channel  thicknesses.  It  also 
generates  the  TACTl  calculational  grid  and  calculates  the  arc  length  distances  between 
nodal  points  in  the  grid  that  are  required  by  TACTl  as  input.  The  geometrical  data  are 
stored  in  a  namelist  data  set  for  use  directly  by  TACTl.  Furthermore,  TACTGRID 
produces  a  scaled  computer  plot  of  each  blade  slice  and  thus  eliminates  the  need  for 
any  intermediate  drafting. 


INTRODUCTION 


The  design  procedure  for  axial-flow  turbine  blading  used  at  the  NASA  Lewis  Re¬ 
search  Center  employs  computer  programs  for  calculating  both  the  flow  field  around  an 
external  blade  shape  and  its  external  heat-transfer  characteristics  (refs.  1  to  3).  For 
blading  with  internal  impingement  cooling,  a  newly  prepared  code  called  TACTl  (ref.  4) 
is  used.  TACTl  is  a  FORTRAN  IV  computer  program  that  calculates  transient  and 
steady-state  temperatures,  pressures,  and  flows  in  a  cooled  turbine  blade  or  vane  with 
an  impingement  insert. 

The  inputs  for  the  computer  programs  of  references  1  to  3  (MERIDL,  TSONIC, 
and  STANS)  are  interconnected;  MERIDL  generates  input  for  TSONIC,  and  TSONIC 
generates  input  for  STANS.  The  geometrical  input  for  MERIDL,  as  well  as  that  for 
TSONIC,  can  be  specified  quite  simply  with  a  few  blade-surface  points.  The  program 
then  uses  spline  curves  passing  through  these  points  to  define  the  blade  contour.  The 
outputs  of  these  programs  cannot  be  directly  input  to  TACTl,  although  TACTl  requires 
output  from  all  three.  The  manual  assembly  of  these  data  for  input  to  TACTl  is  a 
time-consuming  task.  Therefore,  a  computer  program  called  TACTGRID  was  written 
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to  assemble  the  required  data  from  these  codes  and  to  generate  the  internal  geometry 
of  an  impingement-cooled  blade  by  the  use  of  spline  curve  fitting. 

TACTGRID  constructs  the  blade  internal  geometry  from  the  previously  specified 

external  blade-surface  points  and  newly  selected  wall  and  channel  thicknesses .  TACT- 

GBID  also  generates  the  TACTl  calculational  grid  and  calculates  the  arc  length  dis- 

\ 

tances  between  nodal  points  in  the  grid  that  are  required  by  TACTl  as  input.  The  geo¬ 
metrical  data  are  stored  in  a  namelist  data  set  for  use  directly  by  TACTl .  Further¬ 
more,  TACTGRID  produces  a  scaled  computer  plot  of  each  blade  slice  and  thus  elimi¬ 
nates  the  need  for  any  intermediate  drafting.  This  results  in  a  significant  saving  in 
time  and  improvement  in  accuracy . 

The  TACTGRID  program  is  used  at  the  NASA  Lewis  Research  Center  on  an  IBM 
TSS/360-67  computer.  Running  time  for  the  sample  problem,  which  involved  three 
blade  slices,  was  approximately  270  seconds.  The  bulk  of  the  computer  time  is  taken 
up  with  the  computer  plotting . 

The  report  includes  descriptions  of  the  blade  geometric  model,  program  input  and 
output,  and  program  procedure.  The  input  and  output  are  described  in  general  terms 
and  illustrated  by  a  sample  blade  problem . 


DESCRIPTION  OF  TACTl 

The  TACTl  computer  program  was  devised  to  analyze  the  heat-transfer  and 
coolant-flow  distribution  in  an  impingement-cooled  turbine  blade  or  vane .  As  shown  in 
figure  1,  the  coolant  flow  in  this  scheme  enters  the  plenum  through  the  blade  base  and 
flows  sp)anwise  in  the  plenum .  The  coolant  is  discharged  from  the  plenum  in  the  form 
of  jets  from  holes  in  the  plenum  walls.  These  jets  impinge  on  the  inner  surface  of  the 
blade  shell  and  cool  it.  Additional  cooling  is  obtained  as  the  fluid  moves  toward  the 
trailing  edge  through  the  channel  formed  by  the  plenum  wall  and  the  blade  shell. 

In  the  trailing-edge  region  of  the  blade,  the  channel  formed  by  the  blade- shell  in¬ 
ner  surfaces  may  have  pin  fins  for  additional  heat  transfer.  The  coolant  flow  is  ejected 
from  the  blade  through  trailing-edge  holes  or  slots . 

In  order  to  simplify  flow-field  calculations,  the  blade  is  divided  into  slices  (analo¬ 
gous  to  slices  of  bread)  bounded  by  planes  or  surfaces  of  revolution,  as  shown  in  fig¬ 
ure  1 .  These  surfaces  of  revolution  are  placed  a  constant  radial  distance  apart  so  that 
each  slice  has  a  constant  radial  thickness.  A  typical  blade  cross  section  is  depicted  in 
figure  2.  One-dimensional  flow  is  assumed  in  the  flow  channels  shown.  For  each  slice 
an  independent  flow-field  solution  is  performed  that  satisfies  conservation  of  mass, 
momentum,  and  energy  within  that  slice.  In  the  course  of  these  flow-field  calculations, 
heat-transfer  coefficients  are  evaluated  for  two  distinct  regions.  In  the  forward  region, 
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heat  is  transl'erred  by  impingement  cooling  and  forced-convection  cooling  in  the  chan¬ 
nel.  Li  the  trailiag-edge  region,  heat  is  transferred  by  forced-convection  cooling  and, 
when  pin  fins  f.re  specified,  by  the  pin  fins  in  crossflow.  The  blade-shell  temperature 
distribution  is  calculated  from  a  three-dimensional  heat  conduction  model,  which  allows 
for  heat  flow  through  the  wall  and  in  the  spanwise  and  chordwise  directions  as  well. 

Figure  2  shows  the  geometrical  breakdown  of  a  typical  blade  cross  section  into 
calculational  stations.  The  stations  are  odd-numbered  on  the  pressure  side  and  even- 
numbered  on  the  suction  side.  The  inset  in  figure  2  shows  one  station  broken  down  into 
its  five  calculational  nodes,  located  as  follows: 

(a)  A  coating  outer  surface  node 

(b)  A  coating-metal  interface  node 

(c)  A  midmetal  node 

(d)  An  inner  surface  node 

(e)  A  mid-coolant-channel  node 

It  is  evident  that  an  enormous  amount  of  data  must  be  generated  to  completely 
specify  the  blade  geometry  for  use  in  the  TACTl  code.  The  manual  generation  of  these 
data  requires  the  drafting  of  each  blade  cross  section  and  the  measuring  of  coordinates 
and  dimensions .  This  process  is  both  time  consuming  and  subject  to  human  error.  In 
order  to  eliminate  human  error  and  reduce  the  time  required  to  prepare  the  geometri¬ 
cal  input  for  TACTl,  this  process  was  computerized  as  described  in  the  following  sec¬ 
tions  . 


DESCRIPTION  OF  TACTGRID 

The  TACTGRID  program  was  written  to  eliminate  the  manual  drafting  chores  in¬ 
volved  in  designing  a  cooling  configuration  and  running  TACTl .  These  chores  are  now 
performed  by  the  computer,  which  uses  the  ejdernal  blade  profiles  generated  during 
the  aerodynamic  design  process  and  user-specified  distributions  of  coating  thickness, 
metal  thickness,  and  Impingement  crossflow-channel  width.  The  dimensions  and  coor¬ 
dinates  needed  for  TACTl  are  calculated  analytically  with  spline  fit  curves .  There¬ 
fore,  some  restrictions  had  to  be  imposed  on  the  geometry  that  could  be  analyzed  with 
TACTGRID.  The  following  discussion  briefly  describes  the  functions  of  TACTGRID 
and  the  limitations  in  its  use. 

The  internal  geometry  of  each  blade  cross  section  is  generated  by  constructing  in¬ 
ward  normals  to  the  blade  external  surface  at  specified  intervals.  Along  these  nor¬ 
mals,  distances  are  measured  that  correspond  to  the  thicknesses  of  the  coating  (if 
specified),  the  blade  shell,  the  midpoint  of  the  blade  shell,  the  ii.ipingement  crossflow 
channel,  the  midpoint  of  the  channel,  and  the  plenum  wall.  Spline  fit  curves  are  passed 
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through  corresponding  points  on  each  normal.  These  curves  are  called  horizontal  or- 
thogonals  and  variously  represent  the  contours  of  the  interface  between  the  blade  shell 
and  the  coating,  the  outlines  of  the  external  surface  and  the  flow  channels,  and  the  mid¬ 
lines  of  the  blade  shell  and  the  flow  channels.  A  family  of  curves  orthogonal  to  these 
curves,  called  vertical  orthogonals,  is  generated  by  the  predictor-corrector  method  of 
reference  5.  This  method  uses  normals  between  adjacent  horizontal  orthogonals  to 
generate  the  vertical  orthogonals .  The  vertical  orthogonals  represent  the  TACTl  cal¬ 
culation  stations;  and  the  distances  along  them,  between  corresponding  horizontal  or¬ 
thogonals,  are  the  actual  thicknesses  of  coating,  metal,  and  flow  channel  that  are  input 
to  TACTl.  The  distances  along  the  horizontal  orthogonals,  between  adjacent  vertical 
orthogonals,  are  the  actual  chordwise  lengths  input  to  TACTl.  The  procedure  by 
which  TACTGRID  performs  these  tasks  is  described  in  the  Program  Procedure  sec¬ 
tions  that  follow . 

In  regions  of  the  blade  where  the  user-specified  wall  and  channel  thicknesses  have 
no  chordwise  variation,  the  horizontal  orthogonals  constitute  mathematically  parallel 
curves  in  the  plane  of  the  cross  section,  and  the  vertical  orthogonals  are  just  straight 
lines.  The  vertical  orthogonals  wUl  become  curved  lines  in  regions  of  the  blade  where 
the  wall  or  channel  thicknesses  are  not  constant.  TACTGRID  also  calculates,  and 
places  in  the  namelist  input  data  set  for  TACTl,  the  radial  span  for  each  slice. 

Heretofore,  the  finite-difference  grid  for  each  blade  slice,  which  results  from  the 
aforementioned  nodal  arrangement,  had  to  be  generated  by  hand;  a  cross  section  of  a 
(manually  executed)  sample  blade  problem  is  shown  in  figure  3.  TACTGRID  now  con¬ 
structs  this  finite-difference  grid  automatically,  without  any  intermediate  drafting  or 
human  interaction.  TACTGRID  supplies  a  scaled  computer  drawing  of  each  blade  slice 
and  accompanying  grid  as  part  of  its  normal  output  (fig.  4). 


TACTGRID  Limitations 

TACTl  allows  for  mixed  radial-axial  flow  of  the  hot  gas  stream,  so  the  cutting 
planes  (i.e. ,  surfaces  of  revolution)  for  the  slices  do  not  have  to  be  at  constant  radius. 
An  orthogonal  three-dimensional  grid  cannot  be  maintained  as  the  end  walls  become 
greatly  contoured  if  the  grid  lines  are  to  lie  along  the  blade  surface.  Therefore, 
TACTGRID  is  limited  to  constructing  grids  for  blades  with  both  cylindrical  end  walls 
and  cylindrical  surfaces  of  revolution  defining  the  blade  slices . 

TACTGRID  was  developed  specifically  to  automate  the  turbine  blade  design  proce¬ 
dure  at  the  NASA  Lewis  Research  Center.  Since  the  second  step  in  this  procedure  - 
namely,  the  TSONIC  program  (ref.  2)  -  requires  a  specific  t3q>e  of  blade  profile, 
TACTGRID  is  limited  to  generating  grids  for  the  TSONIC  profile.  The  TSONIC  blade 
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shape  and  input  specifications  are  shown  in  figure  5.  The  profile  basically  consists  of 
two  circular  regions  mated  to  two  spline  curves .  TACTGRID  assumes  that  a  TSONIC 
blade  shape  has  been  input  and  uses  the  computer  program  of  reference  6  (TFORM)  to 
transform  that  shape  from  the  TSONIC  coordinate  system  (fig.  5)  to  coordinates  with 
respect  to  true  blade  chord  (figs.  6  to  8).  TACTGRID  neglects  the  trailing-edge  circle 
in  constructing  its  grid. 

In  keeping  with  the  circular  specification  of  the  leading-edge  region  of  the  blade, 
TACTGRID  allows  for  no  wall  thickness  variation  in  the  leading-edge  circular  region 
of  a  slice.  In  this  way,  the  orthogonal  coordinate  lines  that  are  constructed  in  the  nose 
region  can  simply  be  those  of  a  cylindrical  coordinate  system  (fig .  9) .  This  restric¬ 
tion  is  not  present  in  TACTl.  TACTGRID  does,  however,  allow  for  a  completely  gen¬ 
eral  chordwise  variation  of  coating,  metal,  and  impingement-channel  thicknesses  in 
the  spline  regions. 

TACTGRID  does  not  permit  general  variation  of  arc  length  spacing  between  sta¬ 
tions,  although  TACTl  does.  Rather,  the  distribution  is  limited  to  six  regions  with 
uniform  spacing  between  stations  along  the  inner  surface,  within  each  region.  The 
breakdown  of  a  slice  into  six  regions  is  illustrated  in  figure  10.  Equal  increments  of 
inside-wall  arc  length  separate  stations  in  a  given  region  for  the  first  slice .  Succeed¬ 
ing  slices  have  stations  proportionally  alined  with  those  of  the  first  slice.  This  proce¬ 
dure  is  e^qjlained  further  in  the  section  DIFFERENCES  IN  PROGRAM  PROCEDURE 
FOR  SUCCEEDING  SLICES. 

Figure  10  also  shows  that  TACTGRID  will  only  consider  a  straight  trailing-edge 
channel,  whereas  TACTl  does  not  restrict  the  shape  of  the  trailing-edge  channel.  The 
geometry  was  so  limited  because  present  designs  employ  only  drilled  holes  or  rectan¬ 
gular  slots  for  the  trailing-edge  channels.  This  limitation  reduces  fabrication  costs. 

Geometrical  Input  Not  Provided 

TACTl  requires  additional  geometrical  input,  depending  on  the  mode  of  heat  trans¬ 
fer  specified,  that  is  not  provided  by  TACTGRID.  This  input  comprises  the  set  of  vari¬ 
ables  available  to  the  designer  for  optimizing  his  cooling  design.  Since  these  variables 
can  be  obtained  with  no  special  effort,  they  are  not  provided  by  TACTGRID.  These  in¬ 
put  variables  include 

(1)  Impingement-hole  diameter  and  location 

(2)  Pin-fin  diameter  and  location 

(3)  Film-cooling-hole  diameter  and  location 

Figure  11  shows  an  impingement  hole  of  diameter  TDHJ  at  station  ISTA. 
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TACTGRID  does  provide  information  regarding  the  geometrical  co  ’^raints  on  this 
input,  because  it  furnishes  arc  length  distances  along  all  walls .  This  Ci-ables  the  de¬ 
signer  to  prevent  one  station's  row  of  holes  from  overlapping  that  of  an  adjacent  station 
along  the  insert  walls . 


PROGRAM  INPUT 

The  organization  of  a  TACTGRID  input  data  set  for  a  sample  blade  that  was  aero- 
dynamically  designed  for  three  sections  (hub,  mean,  and  tip)  is  illustrated  in  figure 
12(a).  The  input  for  each  blade  slice  contains  data  that  are  needed  solely  for  the 
TACTGRID  program  as  well  as  the  input  data  set  required  by  the  TFORM  program 
(ref.  6),  which  TACTGRID  runs.  The  function  of  the  TFORM  program  in  constructing 
the  TACTl  grid  is  explained  further  in  the  section  PROGRAM  PROCEDURE  FOR  FIRST 
SLICE.  The  TFORM  input  data  set  for  each  slice  follows  the  TACTGRID  input  data  set 
for  that  slice. 


TACTGRID  Input 

An  input  form  for  TACTGRID  is  shown  in  figure  12(b) .  Although  some  input  will  be 
in  real-number  form  (F  format)  and  some  in  integer  form  (I  format),  all  data  will  be 
entered  in  10-column  fields.  In  the  following  description  of  input  variables,  those 
FORTRAN  names  ending  in  T  refer  to  the  suction  side  of  the  blade,  and  those  ending  in 
B  refer  to  the  pressure  side: 

WRU  Width,  in  relative  units,  desired  for  the  scaled  computer  drawing  of  the 

blade  slice.  Each  computer  frame  is  10  relative  units  by  10  relative 
units .  The  plot  is  centered  in  the  frame  so  that  the  left  end  is  at 
(10  -  WRU)/2  relative  units  and  the  right  end  is  at  10  -  (10  -  WRU)/2 
relative  units.  A  WRU  larger  than  9  should  not  be  specified  in  order  to 
maintain  a  0.5-relative-unit  margin  on  either  side  of  the  blade  for  con¬ 
structing  the  labeling  arrows  for  the  nose  and  trailing-edge  stations 
(fig.  4).  TACTGRID  uses  the  IBM  plotting  package  at  the  Lewis  Research 
Center.  TACTGRID  can  be  run  without  the  plotting  by  placing  a  C  in  the 
first  column  of  those  statements  in  the  program  that  call  plotting  subrou¬ 
tines  and  thus  changing  plot  calls  into  comment  cards.  The  names  of  the 
various  plotting  routines  that  TACTGRID  uses  are  listed  in  the  appendix. 
The  format  for  WRU  is  FI 0.5. 
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HRU 


NSLICE 


LU 


NFRAME 


F1T,F2T, 

F1B,F2B 


Height,  in  relative  units,  desired  for  the  scaled  computer  draw¬ 
ing  of  the  blade  slice .  The  drawing  of  the  blade  slice  is  exe¬ 
cuted,  without  axes,  in  a  rectangular  box  of  dimensions  WRU 
relative  units  by  HRU  relative  units.  The  bottom  of  the  slice 
(chord  line)  is  placed  10  -  HRU  relative  units  from  the  bottom  of 
the  computer  frame.  The  scales  on  the  x-  and  y-axes  are  the 
same.  Therefore,  HRU  should  not  be  selected  so  small  as  to 
cause  the  upper  surface  of  the  blade  to  exceed  the  upper  bound¬ 
ary  of  the  frame.  HRU  should  not  be  selected  so  large  as  to 
place  the  chord  line  so  low  in  the  frame  that  there  is  no  room 
left  for  printing  the  legend  and  the  slice  number  below  the  blade 
plot  (fig.  4).  The  format  for  HRU  is  F10.5. 

Number  of  slices  in  the  blade.  The  blade  coordinates  for  each 
slice  are  assumed  by  TACTl  to  be  those  of  a  cross  section  in 
the  middle  of  the  slice.  The  format  for  NSLICE  is  110. 

The  desired  user  units  (output  units)  for  the  problem  are  as 
follows;  If  UU  =  1,  user  ur.its  are  centimeters;  if  UU  =  2,  user 
units  are  inches.  The  format  for  UU  is  110. 

Number  of  frames  desired  for  the  plot.  A  larger  plot  can  be  ob¬ 
tained  by  specifying  NFRAME  larger  than  1.  TACTGRID  will 
multiply  WRU  by  NFRAME  in  determining  the  desired  plot  width 
for  the  blade.  HRU  should  be  adjusted  so  the  blade's  vertical 
dimensions  do  not  exceed  the  borders  of  the  frame.  In  no  in¬ 
stance  can  the  blade  plot  be  higher  than  10  relative  units.  The 
format  for  NFRAME  is  110. 

The  impingement  region  of  the  blade  in  the  spline  portion  is 
divided  into  three  regions  (H,  in,  and  IV)  of  equal  arc-length 
spacing  between  stations  along  the  blade  inside  wall  (fig.  10). 

The  distances  along  the  inside  wall  in  the  impingement  region 
(determined  by  TACTGRID)  are  STWIMP  on  the  suction  side  and 
SBWIMP  on  the  pressure  side.  FIT  is  the  fraction  of  the  dis¬ 
tance  STWIMP  from  the  end  of  region  I  to  the  end  of  region  II  on 
the  suction  side .  F2T  is  the  fraction  of  the  distance  STWIMP 
from  the  end  of  region  I  to  the  end  of  region  HI  on  the  suction 
side.  FIB  is  the  fraction  of  the  distance  SBWIMP  from  the  end 
of  region  I  to  the  end  of  region  II  on  the  pressure  side.  F2B  is 
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NIT,  N2T,  N3T, 
NIB,  N2B,  N3B 


the  fraction  of  the  distance  SBWIMP  from  the  end  of  region  I  to 
the  end  of  region  ni  on  the  pressure  side.  Thus  the  regions  on 
the  suction  and  pressure  sides,  respectively,  are  divided  as 
follows: 

Region  n  =  (F1T*STWIMP),  suction  side 
=  (FIB+SBWIMP),  pressure  side 
Region  III  =  (F2T-F1T)*STWIMP,  suction  side 
=  (F2B-F1B)*SBWIMP,  pressure  side 
Region  IV  =  (1 . -F2T)*STWIMP,  suction  side 
=  (1 . -F2B)*SBWIMP,  pressure  side 
The  format  for  these  variables  is  F10.5. 

Number  of  stations  (vertical  orthogonals)  desired  in  the  circular 
region  (region  I)  on  the  suction  side.  TACTGRID  uses  this  same 
number  of  stations  for  the  pressure  side,  as  TACTGRID  divides 
the  leading-ec^e  circular  region  in  two  and  places  station  1  in 
the  middle.  Station  1  is  common  to  both  the  suction  and  pressure 
sides  and  is  counted  twice.  The  format  for  NVT  is  110. 

Number  of  stations  (vertical  orthogonals)  desired  in  regions  II, 
in,  and  IV  on  the  suction  and  pressure  sides,  respectively. 

These  stations  are  uniformly  spaced  with  respect  to  the  inside¬ 
wall  arc  length  in  each  region.  As  the  last  vertical  orthogonal  in 
the  circular  region  is  not  counted  in  NIT  and  NIB,  the  number  of 
input  stations  equals  the  number  of  arc  length  intervals  between 
stations  in  the  respective  regions .  Since  TACTl  requires  the 
same  number  of  stations  on  the  suction  and  pressure  sides  of  the 
blade,  NIT  +  N2T  +  N3T  must  equal  NIB  +  N2B  +  N3B.  The  for¬ 
mat  for  these  variables  is  110. 

Number  of  stations  desired  in  the  trailing-edge  region  between 
the  impingement  region  and  the  straight  trailing-edge  channel 
region  (region  V  in  fig.  10).  N4B  actually  represents  the  number 
of  stations  (with  equal  inside-wall  arc  length  spacing  between 
them)  on  the  pressure  side  in  region  V.  TACTGRID  determines 
the  locations  of  the  corresponding  suction-side  stations  for  this 
forced- convection,  channel-flow  region.  This  procedure  is 


explained  in  the  section  PROGRA.M  PROCEDURE  FOR  FIRST 
SLICE .  Again,  the  number  of  input  stations  equals  the  number 
of  arc  length  intervals  between  stations  in  this  region.  The  for¬ 
mat  for  N4B  is  110. 

N5B  Number  of  stations  desired  in  the  straight  trailing-edge  channel 

region  (region  VI  in  fig.  10)  on  the  pressure  side.  There  is 
equal  inside-wall  arc  length  spacing  between  these  stations  on 
the  pressure  side.  Again,  the  corresponding  suction-side  sta¬ 
tions  are  determined  by  TACTGRID.  The  number  of  input  sta¬ 
tions  equals  the  number  of  arc  length  intervals  betwe  stations 
in  this  region.  The  format  for  N5B  is  110. 

For  slices  beyond  the  first,  the  preceding  five  lines  of  TACTGRID  input  are 
omitted,  as  they  apply  to  all  the  slices . 

The  wall  and  channel  thicknesses  are  entered  in  namelist  format  for  the  conve¬ 
nience  of  the  user.  The  arrays  for  these  dimensions  are  contained  in  the  TACTGRID 
namelist  name  GEO.  The  distance  along  the  outer  surface  of  the  blade,  for  the  suction- 
and  pressure-side  spline  regions,  is  divided  into  50  equal  increments.  The  elements 
in  the  thickness  arrays  correspond  to  the  51  end  points  of  these  50  arc  length  intervals . 
Thus  a  wall  thickness  distribution  can  be  specified  that  is  based  on  2-percent  incre¬ 
ments  of  spline-region  surface  length  (fig.  9).  The  arrays  that  specify  wall  and  chan¬ 
nel  thicknesses  are  as  follows: 

THKIT,  THKIB  Real-variable  arrays  of  dimension  51  that  specify  the  coating 

thickness  distributions  on  the  suction  and  pressure  sides,  re¬ 
spectively.  The  values  must  be  given  in  the  user-designated 
units.  THKlT(l)  must  equal  THKlB(l). 

THK2T,  THK2B  Real- variable  arrays  of  dimension  51  that  specify  metal  thick¬ 

ness  distributions  on  the  suction  and  pressure  sides,  respec¬ 
tively.  The  values  must  be  given  in  the  user-designated  units. 
THK2T(1)  must  equal  THK2B(1). 

THK3T,THK3B  Real- variable  arrays  of  dimension  51  that  specify  impingement- 

channel  thickness  distributions  on  the  suction  and  pressure  sides, 
respectively .  The  values  must  be  given  in  the  user-designated 
units.  THK3T(1)  must  equal  THK3B(1). 

The  rest  of  the  input  data  is  entered,  once  again,  in  the  usual  real-number  format, 
as  follows: 


THK(l) 


THK(3) 


THKHOL 


THKWAL 

SDUMP 


First  cleinc'iit  in  the  real- variable  ai-ray  TllK,  an  aria^  ol  loui-  ele¬ 
ments.  TlIK(i)  is  the  eoatiiig  thickness  desired  (in  the  user-designated 
units)  at  the  end  ol  the  spline  region  on  the  suction  side.  Alter  cutting 
the  trailing-edge  hole  and  so  determining  i-egion  VI  (lig.  1<';,  T.'.r'r 
GRID  linearly  tapers  the  input  coating  thickjiesses  with  outer  surlace 
arc  length  in  region  VI,  so  that  TI1K1T(51)  THK(l)  (lig.  4).  11  no 

coating  thickness  taper  is  wanted,  the  input  value  lor  TlIK(l)  should  be 
the  same  as  the  input  value  lor  THK1T(51).  The  lormat  loi-  THK(l)  is 
F10.5. 


Third  element  in  the  THK  array.  It  represents  the  coating  thickjiess 
desired  (in  the  user-designated  units)  at  the  end  ol  the  spline  region  on 
the  pressure  side.  The  mput  eoating  thicknesses  are  tapered  linearly 
with  outer  surlace  arc  length  in  region  VI,  so  that  THK1B(51)  -  THK(o( 
(lig.  4).  If  no  coating  thickness  taper  is  wanted,  the  uiput  value  lor 
THK(3)  should  be  the  same  as  the  input  value  lor  THK1B(51).  THK(3) 
and  THK(4)  are  the  metal  wall  thicknesses  on  the  suction  and  pressure 
sides,  respectively,  at  the  end  ol  the  spline  region.  These  values  are 
determined  by  the  program  in  cutting  the  trailing-edge  hole  such  that 
THK(l)  +  THK(2)  =  THK(3)  +  THK(4).  The  lormat  lor  THK(3)  is  Flo. 5. 

Thickness  of  a  straight,  trailing-edge  slot  or  the  diameter  of  a  drilled 
trailing-edge  hole  in  the  user-designat.'d  imits.  If  the  straight,  trail- 
mg-edge  channel  is  intended  to  be  a  drilled  hole  ora  I’ow  of  drilled 
holes,  the  channel  thickness  input  to  TACTl  by  TACTGRID  should  be 
manually  reduced  by  the  user  in  the  following  manner: 

,TllKHOL,^,,CTl - - 

where  JM  is  the  number  of  holes  in  the  slice  and  SPAN  is  the  span  ol 
the  slice.  The  format  for  THKHOL  is  F10.5. 

Thickness  of  the  plenum  wall  in  the  user-designated  units.  The  format 
for  THKWAL  is  F10.5. 


Width  of  the  insert  at  its  end  (fig.  10),  in  the  user-designated  imits. 
The  format  for  SDUMP  is  F10.5. 
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RRO  Radius  of  the  coolant  plenum  outlet  of  the  slice  from  the  axis  of  the 

turbomachine,  in  the  user-designated  units.  The  format  for  RRO  is 
F10.5. 

RRI  Radius  of  the  coolant  plenum  inlet  of  the  slice  from  the  axis  of  the 

turbomachine,  in  the  user-designated  units.  The  format  for  RRI  is 
F10.5. 

Z  Displacement  along  the  TSONIC  z-axis  (fig.  5)  of  the  TSONIC  origin  of 

this  blade  profile  from  the  origin  of  the  profile  in  the  first  slice,  in 
the  user-desigTiated  units.  For  the  first  slice,  Z  is  input  as  0.0.  The 
format  for  Z  is  F10.5. 

THETA  Angular  displacement  of  the  TSONIC  origin  of  this  blade  profile  from 

the  origin  of  the  profile  in  the  first  slice.  THETA  must  be  in  radians. 
For  the  first  slice,  THETA  is  input  as  0.0.  The  format  for  THETA 
is  F10.5. 

In  general,  a  turbine  blade  may  be  designed  with  twist,  taper,  and  lean.  Twist 
refers  to  an  angular  displacement  about  the  stacking  axis  between  successive  blade 
profiles.  Taper  refers  to  a  variation  in  profile  chord  length  along  the  span.  Lean  re¬ 
fers  to  an  angular  displacement  about  the  turbomachine  axis  between  successive  blade 
profiles.  If  the  turbine  blade  design  employs  these  geometric  variations,  the  TSONIC 
origin  (fig.  5)  of  each  blade  profile  may  be  different.  Since  TFORM  sets  up  a  separate 
Cartesian  coordinate  system  for  each  slice,  the  parameters  Z  and  THETA  enable 
TACTGRID  to  transform  the  TFORM  coordinates  for  a  given  slice  back  to  the  TFORM 
coordinate  system  of  the  first  slice,  which  is  taken  as  global. 

TFORM  Input 

An  input  form  for  the  TFORM  input  is  shown  in  figure  12(c).  The  following  de¬ 
scription  of  input  variables  is  excerpted  from  reference  6; 

CHORD  Overall  length  of  blade  in  the  z-direction  (fig.  5) 

STGR  Angular  ©-coordinate  for  center  of  trailing-edge  circle  of  blade  with 

respect  to  center  of  leading-edge  circle  (fig.  5),  rad 

RMI  Radius  of  blade  section  from  axis  of  rotation,  if  RMI  =1,  all  ©-coor¬ 

dinates  are  the  actual  linear  dimension  w. 
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SCALE 


Ratio  of  output  dimensions  to  input  dimensions.  For  example, 
if  input  is  in  feet  and  output  is  desired  in  inches,  SCALE  =  12 
should  be  used. 


DELX 


RI1,RI2 

R01,R02 

BETI1,BETI2 

BET01,BET02 

SPLN01,SPLN02 


MS  PI,  MS  P2 


THSP1,THSP2 


Spacing  of  output  coordinates  in  x-direction  (fig.  8).  DELX 
should  be  chosen  to  be  at  least  CHORD*SCALE/LOO.  DELX 
must  be  given  in  the  output  units  (i.e. ,  if  input  is  in  feet  and 
output  is  in  inches  (SCALE  =  12),  DELX  is  in  inches). 

Leading-edge  radii  of  the  two  blade  surfaces  (fig.  5) 

T railing-edge  radii  of  the  two  blade  surfaces  (fig.  5) 

Angles  (with  respect  to  z-direction)  at  tangent  points  of  lead¬ 
ing-edge  radii  with  the  two  blade  surfaces  (fig.  5),  deg. 

These  must  be  true  angles  in  degrees. 

Angles  (with  respect  to  z-direction)  at  tangent  points  of  trail- 
edge-radii  with  the  two  blade  surfaces 

Number  of  blade  spline  points  given  for  each  surface  as  input 
(maximum,  50) .  These  include  the  first  and  last  points  (dum¬ 
mies)  that  are  tangent  to  the  leading-  and  trailing-edge  radii 
(fig.  5). 

Arrays  of  z-coordinates  of  spline  points  on  the  two  blade  sur¬ 
faces,  measured  from  blade  leading  edges  (fig.  5).  The  first 
and  last  points  in  each  array  must  be  left  blank,  since  these 
values  are  calculated  by  the  program.  If  the  last  point  is  on 
a  new  card,  a  blank  card  must  be  used. 

Arrays  of  ©-coordinates  of  spline  points  corresponding  to 
MSPl  and  MSP2,  rad.  Blanks  must  be  used  in  positions  cor¬ 
responding  to  those  in  MSPl  and  MSP2. 


The  TFORM  input  is  explained  further  in  reference  6.  The  input  data  set  for 
TSONIC  (ref.  2)  is  quite  similar  to  that  for  TFORM  and  provides  all  the  input  needed 
for  TFORM.  The  input  data  set  for  TSONIC  is  generated  by  MEREDL  (ref.  1). 
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PROGRA.M  OUTPUT 


The  output  of  TACTGRID  consists  of  five  parts,  which  are  enumerated  and  de¬ 
scribed  in  the  following  sections . 


TACTGRID  and  TACTl  Overall  Input 

TACTGRID  prints  back  its  input  in  the  labeled  form  illustrated  by  figure  12(b) . 

This  output  is  provided,  for  debugging  purposes,  on  the  first  two  pages  of  the  program 
output  (figs.  13(a)  and  (b)).  The  first  output  page  contains  TACTGRID  overall  input. 
The  second  page  contains  TACTGRID  geometrical  input  for  the  particular  slice.  Also 
contained  on  the  first  output  page  is  the  date  of  the  computer  run  and  the  first  namelist 
input  line  for  TACTl .  The  CHANLS  namelist  name  in  TACTl  pertains  to  the  overall 
job.  For  these  reasons,  the  output  of  figure  13(a)  is  printed  only  here  and  not  repeated 
for  the  various  slices.  The  TACTl  input  is  explained  in  detail  in  reference  4. 


TFORM  Output 

The  TFORM  output  is  provided  on  the  next  three  output  pages  (3  to  5).  The  first 
of  these  pages  (fig.  13(c))  consists  of  a  printout  of  the  TFORM  input  data  set  in  the 
form  of  figure  12(c).  The  second  of  these  pages  (fig.  13(d))  contains  blade  data  at  the 
input  spline  points.  This  information  is  used  for  debugging  because  wildly  varying  sec¬ 
ond  derivatives  at  the  input  spline  points  will  indicate  an  input  error.  The  third  of 
these  pages  (fig.  13(e))  gives  the  transformed  blade  coordinates.  The  first  line  gives 
the  output  number  of  x-coordinates  and  the  orientation  angle  (p  (fig.  8).  The  next  lines 
present  a  tabulation  of  the  x-  and  y-coordinates  for  the  upper  (suction  side)  and  lower 
(pressure  side)  surfaces  (fig.  8).  Any  TFORM  error  messages  will  also  be  printed  out. 


TACTl  Input 

On  the  sixth  output  page  (fig.  13(f))  begins  the  TACTl  input  pertaining  to  the  indi¬ 
vidual  blade  slice.  The  TACTl  namelist  names  for  which  TACTGRID  provides  input 
are  CHANLS,  CONTRL,  PROPS,  and  GEO.  Within  each  namelist  name,  TACTGIJD 
only  provides  the  input  for  those  variables  that  are  geometry  related.  AU  variables 
not  so  related  are  omitted  from  the  namelist  name  printout .  The  variables  provided, 
within  their  respective  namelist  names,  are  as  follows: 
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Namelist  /CHANLS/  NS  LICE,  NSTA,  lUNITS 
Namelist  /CONTRL/  NFWD,  NGEO 

Namelist  /PROPS/  SPAN,  ADUMP,  DHYD,  APLEN,  RI,  RO 
Namelist  /GEO/  ISTA,  BTB,  THK,  TDLX 

The  value  for  ADUMP  that  the  program  calculates  is  the  cross-sectional  area  of 
the  plenum  opening  at  its  end,  for  the  "cut-off"  end  shape  assumed  by  TACTGRID 
(fig.  4).  For  designs  employing  other  plenum  end  geometries,  this  input  to  TACTl 
should  be  modified  accordingly . 


Nodal  and  Plenum  Contour  Coordinates 

Representative  blade  coordinate  output  is  shown  in  figure  13(g) .  TACTGRID  calcu¬ 
lates  the  arc  length  distance  along  each  horizontal  contour  line  (horizontal  orthogonal) 
from  an  initial  radial  line  that  is  the  TSONIC  z-axis  (fig.  14).  Distances  clockwise  of 
this  line  are  positive  and  those  counterclockwise  are  negative.  TACTGRID  also  calcu¬ 
lates  the  arc  length  distances  between  nodes  along  the  vertical  orthogonals  joining  them. 
Horizontal  and  vertical  arc  length  coordinates  as  well  as  x,  y-coordinates  (the  x,y- 
coordinates  of  the  TFORM  coordinate  system  set  up  for  the  first  slice,  fig.  7)  are 
given  for  every  grid  node  and  every  intersection  point  of  a  station  vertical  orthogonal 
with  the  insert  outer  and  inner  contours.  In  the  arrays  of  output  coordinates,  the  first 
index,  I,  indicates  the  vertical  orthogonal  line  and  the  second  index,  J,  indicates  the 
horizontal  orthogonal  line.  All  output  coordinates  are  in  the  user-designated  units. 

The  horizontal  orthogonals  are  denoted  as  follows: 

J  =  1  represents  the  blade  outer  surface. 

J  =  2  represents  the  coating-metal  interface  line. 

J  =  3  represents  the  midmetal  line. 

J  =  4  represents  the  blade  inner  surface. 

J  =  5  represents  the  impingement-channel  midchannel  line. 

J  =  6  represents  the  insert  outer  surface. 

J  =  7  represents  the  insert  inner  surface. 
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The  output  coordinates  are  as  follows  (fig.  14): 


XCT(I,J),XCB(I,J) 


YCT(I,J),YCB(I,J) 


SCT(I.J),SCB(I.J) 


NCT(I,J),NCB(I,J) 


XTT(I,J),XBB(I,J) 


YTT(I,J),YBB(I,J) 


Arrays  of  x-coordinates  of  grid  points  in  the  circular  region 
on  the  suction  and  pressure  sides,  respectively.  The 
TSONIC  z-axis  and  the  borders  of  the  circular  region  are 
delineated  by  radial  lines  on  all  output  blade  plots  (figs.  14 
and  4).  1  =  1  corresponds  to  station  1  for  both  sides.  1  =  4 
refers  to  the  fourth  station  in  the  circular  region  (station  1 
included)  on  the  respective  side.  For  instance,  XCT(4,  3) 
is  the  x-coordinate  of  the  midmetal  node  in  station  6  of  fig¬ 
ure  4(a).  The  dimension  of  these  arrays  is  10  x  7. 

Arrays  of  y-coordinates  of  grid  points  in  the  circular  region 
on  the  suction  and  pressure  sides,  respectively.  The  dimen¬ 
sion  of  these  arrays  is  10  x  7. 

Arrays  of  horizontal  arc  length  coordinates  in  the  circular 
region  on  the  suction  and  pressure  sides,  respectively.  The 
dimension  of  these  arrays  is  10  x  7 . 

Arrays  of  vertical  arc  length  coordinates  in  the  circular  re¬ 
gion  on  the  suction  and  pressure  sides,  respectively.  Here 
the  J  index  refers  to  the  distance  between  the  outer  coating 
node  and  the  J  +  1  node.  For  instance,  NCT(4, 3)  is  the  dis¬ 
tance  along  the  vertical  orthogonal  of  station  6  from  the 
outer  coating  node  to  the  metal  inner  surface  node  of  figure 
4(a) .  The  dimension  of  these  arrays  is  10  x  7 . 

Arrays  of  x-coordinates  of  grid  points  in  the  spline  region 
on  the  suction  and  pressure  sides,  respectively.  1  =  1  cor¬ 
responds  to  the  first  station  in  the  spline  region.  1=4  re¬ 
fers  to  the  fourth  station  in  the  spline  region  on  the  respec¬ 
tive  side.  For  instance,  XTT(4,  3)  is  the  x-coordinate  of  the 
midmetal  node  in  station  14  of  figure  4(a) .  The  dimension  of 
these  arrays  is  50  x  7 . 

Arrays  of  y-coordinates  of  grid  points  in  the  spline  region 
on  the  suction  and  pressure  sides,  respectively.  The  dimen¬ 
sion  of  these  arrays  is  50  x  7 . 
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STT(I,  J),SBB(I,  J)  Arrays  of  horizontal  arc  length  coordinates  in  the  spline  re¬ 

gion  on  the  suction  and  pressure  sides,  respectively.  The 
dimension  of  these  arrays  is  50  x  7. 

NTT(I,  J),  NBB(I,  J)  Arrays  of  vertical  arc  length  coordinates  in  the  spline  re¬ 

gion  on  the  suction  and  pressure  sides,  respectively.  As  in 
the  circular  region,  the  J  index  refers  to  the  distance  be¬ 
tween  the  outer  coating  node  and  the  J  +  1  node.  For  in¬ 
stance,  NTT(4, 3)  is  the  distance  along  the  vertical  orthogo¬ 
nal  of  station  14  from  the  outer  coating  node  to  the  metal  in¬ 
ner  surface  node  of  figure  4(a) .  The  dimension  of  these  ar¬ 
rays  is  50  X  7 . 

Figure  14  illustrates  the  output  coordinate  designations  for  a  few  stations .  The 
TSONIC  computer  program  (ref.  2)  provides  the  boundary  condition  (inviscid  tangential 
velocity  distribution)  for  the  STAN5  boundary- layer  program  (ref.  3)  as  a  function  of 
surface  length  from  the  TSONIC  origin.  Therefore,  the  heat-transfer-coefficient  dis¬ 
tribution,  or  heat-flux  distribution,  which  is  output  by  STAN5,  will  likewise  be  in 
terms  of  surface  length  measured  from  the  TSONIC  origin.  Since  the  TACTl  program 
requires  that  the  heat-transfer  boundary  condition  be  specified  in  terms  of  surface 
length  from  station  1,  the  quantity  SCT(1, 1)  (or  equivalently  SCB(1, 1))  must  be  sub¬ 
tracted  from  each  arc  length  entry  when  STANS  output  is  used  as  input  for  TACTl . 

Scaled  Computer  Plot  of  Blade  Slice 

TACTGRED  produces  a  scaled  computer  drawing  of  every  blade  section  used  for  a 
TACTl  run.  The  designer  can  select  all  the  input  variables  shown  in  figure  12(b)  with¬ 
out  first  making  a  blade  layout  on  a  drafting  board.  By  a  visual  inspection  of  the  com¬ 
puter  drawing,  the  designer  can  readily  detect  gross  errors  in  the  geometrical  dimen¬ 
sions  selected.  The  labeling  on  each  drawing  includes  the  profile  designation  (title 
from  the  TFORM  input  data  set)  in  the  top  left  comer,  the  date  of  the  run  in  the  top 
right  corner,  the  scale  and  legend  in  millimeters  and  inches  centered  below  the  profile, 
and  the  slice  number  centered  on  the  bottom  of  the  frame  (fig.  4).  Every  calculation 
station  is  indicated  with  an  arrow  and  the  appropriate  station  number .  Since  the  draw¬ 
ing  is  scaled,  the  designer  can  plan  and  draw  in  his  design  variables  (impingement 
holes,  pin  fins,  film-cooling  holes)  right  on  the  plot,  with  drafting  implements  and  an 
engineering  scale. 
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PROGRAM  PROCEDURE  FOR  FIRST  SLICE 


For  each  slice,  TACTGRID  first  runs  TFORM  to  change  the  TSONIC  external 
blade  coordinates  from  z, 6 -coordinates  (fig.  5)  to  x, y-coordinates,  where  the  x-axis 
is  tangent  to  the  blade  lower  surface  (fig.  7).  TFORM  obtains  y-coordinates  on  the 
blade  upper  and  lower  surfaces  for  equal  increments  in  x  (fig.  8).  TACTGRID  then 
plots  these  x,y-points  so  that,  if  the  program  is  being  run  interactively  with  an  on-line 
cathode  ray  tube  oscilloscope  (CRT),  the  analyst  can  detect  an  input  error  visually. 

Next,  TACTGRID  calculates  the  x,  y-coordinates  of  the  circle-spline  juncture 
points,  which  were  determined  in  z,0 -coordinates  by  TFORM  (fig.  5),  in  order  to  re¬ 
solve  the  profile  into  a  spline  region  and  a  circular  region  (fig.  10).  A  different  pro¬ 
cedure  is  used  for  generating  the  grid  in  each  region. 

TACTGRID  generates  the  grid  in  the  spline  region  by  first  calculating  the  spline- 
region  outer  surface  length  on  the  suction  and  pressure  sides.  As  mentioned  earlier, 
TACTGRID  divides  these  spline- region  surface  lengths  into  50  equal  increm'  .its  and 
determines  the  51  end  points  of  these  50  increments . 

TACTGRID  then  calculates  all  the  nodal  points  in  the  blade  shell  for  the  circular 
region  and  plots  this  region's  horizontal  and  vertical  orthogonals.  The  orthogonal  grid 
in  the  circular  region  is  just  a  polar  coordinate  system,  with  the  vertical  orthogonals 
being  radial  lines  and  the  horizontal  orthogonals  being  sections  of  circles.  The  pro¬ 
gram  determines  the  positions  of  the  nodal  points  in  the  circular  region  by  the  number 
of  vertical  orthogonals  (stations)  input  because  stations  are  placed  at  equal  angular  in¬ 
crements  (fig.  9(b)). 

TACTGRID  then  defines  the  blade  internal  shape  in  the  spline  region  (fig.  9(a))  by 
generating  normals  to  the  outer  contour  in  this  region  at  the  51  end  points  previously 
determined.  The  coating-metal  interface  line,  on  the  suction  side,  is  formed  by  points 
at  distances  THKIT(I)  along  the  51  normals.  The  midmetal  line  is  formed  by  points  at 
distances  THKIT(I)  +  (THK2T(I)/2)  along  these  normals,  and  the  metal  inside  wall  is 
formed  by  points  at  distances  THKIT(I)  +  THK2T(I)  along  the  normals  (fig.  9(b)).  The 
same  procedure  is  followed  on  the  pressure  side.  Generally,  the  suction-  and  pres¬ 
sure-side  inner  walls  will  intersect  (fig.  9(a)). 

At  this  point,  TACTGRID  prepares  to  cut  a  straight,  trailing-edge  channel  of  thick¬ 
ness  THKHOL  (input)  through  the  generated  blade  shell  (fig.  9(a))  by  an  iterative  proce¬ 
dure.  Because  this  straight  channel  is  to  represent  a  drilled  hole  or  rectangular  slot, 
cut  inwards  from  the  trailing  edge,  the  program  employs  the  same  reasoning  as  a  ma¬ 
chinist  might  apply  in  determining  the  angle  of  the  cut.  Logical  criteria  for  this  cut  are 

(1)  The  location  of  the  channel  should  be  such  that  the  suction-  and  pressure-side 
wall  thicknesses  at  the  trailing-edge  exit  are  equal  (i.e. ,  lines  EB  and  FD  in  figs.  15 
and  16  should  be  equal) . 
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(2)  The  sides  of  the  channel  should  intersect  the  suction-  and  pressure-side  inner 
walls  at  points  that  are  equidistant  from  the  axis  of  the  channel  (points  A  and  I  in  fig. 
16(b)). 

The  program  determines  the  orientation  of  the  channel  by  a  trial-and-error  method 
consisting  of  two  levels  of  iteration.  The  procedure  is  to  pick  a  point  on  the  suction- 
side  inner  surface  as  a  first  guess  for  the  intersection  point  of  the  channel  top  surface 
with  the  suction-side  inner  wall.  After  adjusting  the  angle  of  the  channel  axis  so  that 
the  normals  to  the  external  surface  at  the  trailing  edge  are  cut  at  equal  lengths,  the 
program  checks  to  see  if  the  channel  intersection  with  the  pressure-side  inner  surface 
lies  directly  across  the  channel  from  the  suction-side  intersection  point.  If  it  does  not, 
a  new  suction-side  inner  surface  point  is  selected,  and  so  on.  The  details  of  this  pro¬ 
cedure  are  as  follows  (figs.  15  and  16): 

(1)  Outer  loop;  Pick  point  A  on  the  suction-side  inner  surface  and  point  B  on  the 
suction-side,  trailing-edge  normal  ^  (fig.  15(a))  and  enter  the  inner  loop. 

(2)  Inner  loop:  Construct  a  normal  of  length  THKHOL  to  line  AB  at  point  B,  and 
locate  point  C  at  the  terminus  of  this  normal  (fig.  15(a)). 

(3)  Inner  loop:  Construct  a  line  parallel  to  AB  from  point  C,  and  determine  the 

intersection  point,  D,  of  this  line  with  the  pressure-side,  trailing-edge  normal  FH 
(fig.  15(a)).  _  _  _ 

(4)  Inner  loop:  Calculate  lengths  |EB|  and  |FD|  .  If  |]^|  |FD|  ,  move  point  B 
along  EG  and  return  to  step  (2).  If  |EB|  =  1FD|  (i.e.,  |  lEB|  -  |FDl  |/thKHOL  < 
0.01),  go  to  step  (5)  (return  to  outer  loop).  Figures  15(a)  and  (b)  illustrate  the  first 
and  fifth  iterations,  respectively,  of  the  inner-loop  procedure  for  the  first  outer-loop 
iteration. 

(5)  Outer  loop:  Extend  line  CD  (parallel  to  AB)  to  intersect  the  pressure-side 
inner  wall,  and  determine  the  intersection  point  I  (fig.  16(a)). 

(6)  Outer  loop:  Construct  a  normal  of  length  THKHOL  to  line  AB  at  point  A,  and 
locate  point  J  at  the  terminus  of  this  normal  (fig.  16(a)).  I  and  J  must  coincide  to  meet 
the  cutting  criteria. 

(7)  Outer  loop:  By  construction,  I  and  J  are  colinear  (fig.  16(a)).  If  the  length 
IIJl  0,  move  point  A  along  the  suction-side  inner  surface  and  return  to  step  (2).  If 
|IJ|  =  0  (i.e. ,  llJ|/THKHOL  <  0.01),  exit  from  the  loop  because  the  solution  has  been 
found.  Figures  16(a)  and  (b)  show  the  first  and  17th  iterations,  respectively,  of  the 
outer- loop  procedure  for  cutting  the  trailing-edge  hole. 

At  this  point,  TACTGRID  tapers  the  coating  thickness  linearly,  with  outer  surface 
arc  length,  from  the  value  found  at  the  beginning  of  region  VI  (fig.  10)  to  the  value 
specified  by  the  designer  for  the  trailing  edge  of  the  blade  (THK(l)  on  the  suction  side 
and  THK(3)  on  the  pressure  side) .  This  procedure  redetermines  the  second  horizontal 
contour  (coating-metal  interface  line). 
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The  program  then  recomputes  the  coordinates  of  the  inner  wall  and  the  midmetal 
line  because  of  the  material  removed  in  cutting  the  trailing-edge  channel. 

TACTGRID  now  constructs  the  three  horizontal  contours  needed  for  the  impinge¬ 
ment  crossflow  channel  and  insert  walls  (fig.  4)  by  generating  normals  to  the  blade- 
shell  iimer  surface  at  the  end  points  of  the  outer  surface  normals.  The  first  horizontal 
contour  (impingement-channel  mean  streamline)  is  formed  by  points  at  distances 
THK3T(I)/2  (suction  side)  and  THK3B(I)/2  (pressure  side)  along  the  blade-shell  inner 
surface  normals.  The  second  horizontal  contour  (insert  outer  surface)  is  formed  by 
points  at  distances  THK3T(I)  (suction  side)  and  THK3B(I)  (pressure  side)  along  the  in¬ 
ner  surface  normals.  The  third  contour  (insert  inner  surface)  is  formed  by  points  at 
distances  THK3T(I)  +  THKWAL  (suction  side)  and  THK3B(I)  +  THKWAL  (pressure  side) 
along  the  normals  (fig.  9(b)). 

The  program  places  calculational  stations  at  the  desired  arc  length  positions  desig¬ 
nated  by  the  input  parameters  (FIT,  F2T,  FIB,  F2B,  NIT,  N2T,  N3T,  NIB,  N2B, 

N3B,  N4B,  and  N5B;  see  section  TACTGRID  Input)  by  calculating  the  spline  arc  length 
of  the  inside  wall  (for  the  suction  and  pressure  sides).  At  this  point,  TACTGRID  plots 
the  inside-wall  contours  to  provide  a  quick  screening  check  for  the  designer  who  is 
running  conversationally  with  a  CRT.  After  a  visual  inspection,  the  user  may  decide 
to  terminate  the  program  and  rerun  with  a  different  set  of  input  parameters. 

TACTGRID  now  prepares  to  construct  the  orthogonal  TACTl  grid.  First,  the 
chordwise  extent  of  the  insert  is  determined  by  an  iterative  procedure  (the  details  of 
which  are  described  later).  This  procedure  employs  the  construction  of  vertical  or¬ 
thogonal  links  by  the  method  of  reference  5.  Once  the  impingement  region  has  been 
thus  delimited,  the  nodal  points  on  the  inside  wall  are  located  as  described  in  the  sec¬ 
tion  TACTGRID  Input.  From  these  points  on  the  inner  surface,  vertical  orthogonals 
are  generated  in  the  blade  shell  that  cut  at  right  angles  to  all  the  horizontal  contours 
previously  generated.  The  horizontal  contours  thus  constitute,  in  effect,  horizontal 
orthogonals  for  the  grid.  As  before,  the  method  of  generating  the  vertical  orthogonals 
is  that  employed  in  reference  5. 

The  construction  of  a  vertical  orthogonal  link  by  the  aforementioned  method  is  de¬ 
picted  in  figure  17  (fig.  5  of  ref.  5).  In  this  figure,  a  vertical  "orthogonal"  link  is  gen¬ 
erated  between  two  horizontal  orthogonals  in  an  analogous  r,  z  rectangular  coordinate 
system.  Beginning  at  nodal  point  (REFZ,REFR)  on  the  lower  horizontal  orthogonal,  a 
normal  is  constructed  (line  1  in  fig.  17),  and  the  intersection  (ZNL,  RNL)  is  calculated 
on  the  upper  horizontal  curve.  From  the  slope  SLU  at  this  intersection  point,  line  2  in 
figure  17  is  constructed  in  such  a  way  that  it  is  perpendicular  to  the  SLU  line  and  also 
passes  through  (REFZ,  REFR)  on  the  lower  horizontal  orthogonal.  After  the  coordi¬ 
nates  (ZNU,  RNU)  of  the  intersection  of  line  2  with  the  SLU  line  are  calculated,  the 
final  vertical  orthogonal  link  is  determined  as  the  line  that  intersects  the  upper  curve 
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at  an  abscissa  ZOM  =  (ZNL  +  ZNU)/2.  This  is  the  point  (ZOM,  ROM)  in  figure  17.  De¬ 
tails  of  the  procedure  are  given  in  reference  5.  TACTGRID  uses  this  same  method  to 
construct  approximate  "streamlines”  and  streamline  normals  in  the  impingement 
cixjssflow-channel  region  (regions  II,  III,  and  IV  of  fig.  10)  and  in  the  trailing-edge 
flow  region  (regions  V  and  VI  of  fig.  10). 

TACTGRID  locates  the  end  of  the  insert  within  the  blade  shell,  and  thus  delimits 
the  impingement  region,  by  the  following  procedure  (fig.  18):  Beginning  with  the  point 
on  the  pressure-side  insert  outer  wall  determined  by  the  25th  impingement-channel 
normal,  TACTGRID  constructs  vertical  "orthogonal"  links,  as  previously  described, 
between  the  pressure-side  insert  outer  wall  and  the  suction-side  insert  outer  wall. 

This  point  on  the  pressure-side  insert  outer  wall  is  moved  back  and  forth  in  a  bisection 
root -finding  procedure  until  the  length  of  the  generated  vertical  orthogonal  link  equals 
the  specified  end  width  of  the  insert,  SDUMP.  The  first  four  iterations  of  this  proce¬ 
dure  (lines  1  to  4  in  fig.  18)  as  well  as  the  last  iteration  (21st  iteration  -  line  labeled 
"SDUMP")  are  shown  in  figure  18. 

The  practice  of  constructing  so-called  "approximate  streamlines"  in  the  coolant 
flow  regions  is  motivated  by  a  desire  to  minimize  the  errors  inherent  in  accepting  a 
one-dimensional  flow  approximation.  By  approximating  the  streamlines,  TACTGRID 
is  able  to  construct  cross  sections  that  are  more  nearly  normal  to  the  flow.  The  outer 
streamlines  in  the  impingement  crossflow  channel  are  taken  as  the  channel  boundaries, 
and  a  mean  streamline  is  determined  by  the  spline  curve  passing  through  the  half- 
length  points  along  the  inner  surface  normals,  as  previously  explained.  In  the  trailing- 
edge  region,  vertical  "orthogonal"  links  are  constructed  between  the  pressure-side  in¬ 
ner  wall  and  the  suction-side  inner  wall,  from  the  previously  determined  station  posi¬ 
tions  on  the  pressure-side  inner  wall.  Nine  points  are  taken  equidistant  along  these 
vertical  links.  The  horizontal  spline  curves  that  pass  through  these  points  comprise 
the  approximate  streamlines  (fig.  19).  TACTGRID  then  constructs  eight  vertical  or¬ 
thogonal  links  across  the  approximate  streamlines  to  determine  the  channel  cross- 
sectional  flow  areas.  Figure  20  shows  the  resulting  vertical  orthogonals  and  the 
trailing-edge-channel  mean  streamline  (dotted  line),  which  is  constructed  by  passing 
one  spline  curve  between  each  pair  of  vertical  orthogonals  at  their  midpoints .  Each 
horizontal  spline  curve  between  two  vertical  orthogonals  is  normal  to  the  vertical  or¬ 
thogonals  at  its  end  points.  Figure  20  shows  that  in  TACTGRID  the  array  elements 
NTRLDG(1, 9)  play  the  roles  of  THK(3)  in  TACTl  and  that  the  differences  STRLDG(I)  - 
STRLDG(I-l)  are  used  to  calculate  TDLX(5)  for  TACTl.  The  points  where  the  flow- 
passage  vertical  orthogonals  intersect  the  suction-side  inner  surface  in  regions  V  and 
VI  determine  the  station  positions  on  the  suction  side  in  these  regions. 

If  the  coating,  metal,  and  impingement-channel  thicknesses  have  no  chordwise 
variation,  the  horizontal  contours  previously  mentioned  will  constitute  mathematically 
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parallel  curves  in  the  plane  of  the  cross  section  (except  for  the  wall  portions  in  region 
VI),  and  the  vertical  orthogonals  will  be  straight  lines.  This  is  a  so-called  "normal" 
coordinate  system  that  is  commonly  employed  in  two-dimensional  boundary  layer  cal¬ 
culations.  However,  if  chordwise  variation  is  specified,  and  also  because  of  the  auto¬ 
matic  wall  thickness  tapering  in  region  VI,  a  situation  will  be  encountered  wherein  the 
vertical  orthogonals  are  curved.  The  midmetal  line  and  impingement-channel  mean 
streamline,  when  constructed  as  has  been  e^lained,  will  not  cut  these  vertical  orthog¬ 
onals  precisely  in  two.  Nonetheless,  TACTl  requires  that  the  vertical  orthogonals  be 
divided  into  two  equal  segments.  Therefore,  TACTGRID  has  built  into  it  a  procedure 
that  will  automatically  redefine  the  midmetal  horizontal  orthogonals  and  reconstruct  the 
vertical  orthogonals,  whenever  the  vertical  orthogonals  are  not  split  precisely  in  two. 
This  is  an  iterative  procedure  that  alters  the  horizontal  orthogonals  and  then  recon¬ 
structs  the  vertical  orthogonals,  over  and  over  again,  until  the  vertical  orthogonals  are 
equally  divided.  Experience  with  the  program  has  shown  that  when  no  chordwise  thicx- 
ness  variation  is  specified,  the  midmetal  and  mean- stream  line  contours  will  split  the 
vertical  orthogonals  closely  enough  in  region  VI  that  the  iterative  procedure  is  never 
triggered. 

TACTGRID  calculates  the  arc  lengths  between  nodal  points  along  the  horizontal  and 
vertical  orthogonals  after  the  vertical  orthogonals  have  been  constructed.  All  the  hori¬ 
zontal  contours,  as  well  as  the  vertical  orthogonals  in  the  trailing-edge  flow  region, 
are  represented  mathematically  by  spline  curves  which  are  obtained  by  using  the  curve¬ 
fitting  programs  of  references  1  and  2.  TACTGRID  calculates  the  perimeter  of  the  in¬ 
side  wall  of  the  insert,  SPLNUM,  as  well  as  the  cross-sectional  area  of  the  interior  of 
the  plenum,  APLNUM.  The  cross-sectional  area  is  calculated  by  using  the  formula  for 
the  area  underneath  a  spline  curve  from  subroutine  INTGRL  of  reference  7.  The  hy¬ 
draulic  diameter  of  the  insert,  which  is  provided  to  TACTl  along  with  the  cross-sec¬ 
tional  area,  is  calculated  as  DHYDRP  =  4+APLNUM/SPLNUM . 

DIFFERENCES  IN  PROGRAM  PROCEDURE  FOR  SUCCEEDING  SLICES 

TACTGRID  employs  a  number  of  modifications  to  this  procedure  when  the  slice 
number,  ISLICE,  is  greater  than  1. 

In  the  circular  region,  if  ISLICE  is  greater  th.'^n  1,  TACTGRID  skips  the  station 
position  calculations  based  on  the  number  of  vertical  orthogonals  selected  for  this  re¬ 
gion  and  instead  determines  station  positions  by  placing  the  stations  at  the  same  frac¬ 
tions  of  external- surface  arc  length  that  they  occupied  in  the  first  slice.  This  is  to 
prevent  severe  misalinement  of  stations  between  slices,  which  would  invalidate  the 
orthogonality  assumption  implicit  in  TACTl .  If  some  of  the  calculation  stations  that 
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were  on  the  circle  in  the  first  slice  fall  beyond  the  circular  region  in  the  new  slice,  the 
vertical  orthogonals  for  these  stations  are  constructed  like  those  for  spline-region  sta¬ 
tions  . 

In  the  spline  region  also,  if  ISLICE  is  greater  than  1,  station  positions  are  deter¬ 
mined  by  placing  the  stations  at  the  same  fractions  of  external-surface  arc  length  that 
they  occupied  in  the  first  slice.  The  station  vertical  orthogonals  are  now  generated  in 
the  blade  shell  from  these  points  on  the  outside  wall,  rather  than  from  points  on  the  in¬ 
side  wall  as  was  done  for  the  first  slice.  The  station  positions  on  the  suction-side  in¬ 
ner  surface  for  regions  V  and  VI  are  determined  as  before,  however;  and  the  vertical 
orthogonals  are  generated  from  the  inside  wall  in  regions  V  and  VI  for  the  suction  side. 
If  some  of  the  spline-region  calculation  stations  of  the  first  slice  fall  in  the  circular 
region  in  the  new  slice,  the  vertical  orthogonals  for  these  stations  are  constructed  like 
those  for  circular-region  stations. 

The  position  of  the  plenum  end  is  located  as  before,  but  now  the  end  will  no  longer 
line  up  with  calculational  stations  (figs.  4(b)  and  (c)).  The  stations  marking  the  end  of 
the  impingement  region  will  be  considered  as  the  last  stations  on  the  suction  and  pres¬ 
sure  sides  before  the  plenum  end.  The  output  coordinates  of  the  plenum  end  are  print¬ 
ed  out  even  for  slices  beyond  the  first.  However,  for  ISLICE  greater  than  1,  the  first 
index  of  the  output  coordinates  for  these  four  end  points  is  NPLNUM  +  1  rather  than 
NPLNUM,  where  NPLNUM  is  the  sum  of  the  number  of  stations  in  regions  n,  HI,  and 
IV  for  the  particular  slice.  For  instance,  the  horizontal  arc  length  coordinate  of  the 
plenum  end  point  on  the  suction-side  insert  outer  surface,  in  figure  4(b),  will  be  listed 
in  the  output  as  STT(17,6),  because  there  are  16  impingement- region  stations  in  the 
spline  region  on  the  suction  side. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  August  23,  1979, 
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APPENDIX  -  PROGRAM  STRUCTURE  AND  FUNCTION 


The  TACTGRID  program  consists  of  a  main  program,  GRIDIN,  and  seven  primary 
subroutines  (SGRID2,  SGRIDa,  SGRID4,  SGRID5,  SGRID6,  SGRID7,  andSGRIDB). 
These  subroutines,  in  turn,  call  a  number  of  special-purpose  subroutines  (secondary 
subroutines).  The  calling  relations  between  GRIDIN  and  the  various  subroutines  are 
shown  in  figure  21.  The  vertical  positions  of  the  subroutines  in  figure  21  reflect  their 
order  of  occurrence  in  the  calling  sequence.  They  are  described  in  this  order  here. 

As  shown  in  figure  21,  the  flow  of  the  TACTGRID  program  is  divided  among  the 
main  program  GRIDIN  and  subroutines  SGRID2  to  SGRID8  in  a  linear  progression  of 
the  grid  construction.  This  was  done  merely  to  accommodate  the  compilation  storage 
limitations  of  the  IBM  TSS/360-67  computer  at  the  Lewis  Research  Center,  on  which 
the  TACTGRID  program  is  mounted.  The  functions  of  GRIDIN  and  each  of  the  subrou¬ 
tines  in  TACTGRID  are  described  in  this  appendix.  The  word  descriptions  of  subrou¬ 
tines  obtained  from  other  sources  are  taken  from  the  references  cited. 


Main  Program  GRIDIN 

The  main  program  GRIDIN  reads  in,  and  prints  back  out,  all  TACTGRID  input.  It 
calls  subroutine  TIME  to  obtain  the  date  and  time  of  the  run.  GRIDIN  contains  those 
plot  labels  that  are  common  to  all  the  blade  slice  plots.  It  initializes  the  counters  for 
the  number  of  stations  in  the  circular  region,  the  impingement  spline  region,  and  the 
trailing-edge  spline  region  (for  the  suction  and  pressure  sides)  as  the  values  for  the 
first  slice.  These  values  will  be  altered  by  subroutines  GRID2N  and  GRID4N  for  suc¬ 
ceeding  slices. 

GRIDIN  calls  subroutine  TFORM  to  transform  the  external  blade  contour  descrip¬ 
tion  from  (z,0)-coordinates  to  Cartesian  (x,y)  coordinates.  It  converts  all  geometri¬ 
cal  quantities  that  were  not  previously  converted  by  TFORM  from  (z,0)-coordinates  to 
(x,y)-coordinates.  It  also  converts  all  those  variables  still  in  the  input  units  to  the 
user-designated  units.  It  retains  the  parameters  of  the  first  slice's  coordinate  system 
for  converting  later  slice  output  coordinates  to  the  first  slice' s  coordinate  system, 
which  is  taken  as  global. 

GRIDIN  sets  up  the  data  set  BLADE .  PLOTS  to  store  the  plots  of  the  blade  slices. 
It  sets  up  the  identification  number,  the  imaginary  axes,  and  the  scale  for  each  slice 
plot.  GRIDIN  plots  the  (x,y)-coordinates  returned  by  TFORM.  It  also  plots  the  two 
leading-ec^e  circle  centers  (the  circular  regions  on  the  suction  and  pressure  sides 
need  not  have  a  common  origin;  ref,  4),  the  origin  of  the  TSONIC  coordinate  system 
(same  as  TFORM  origin),  and  the  line  joining  these  points.  This  line  is  the  z-axis  of 
the  TSONIC  coordinate  system. 
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GRIDIN  calculates  the  angular  extent  of  the  suction-  and  pressure-side  circular 
regions  (from  the  TSONIC  axis)  and  the  slope  of  the  profile  at  the  four  circle-spline 
juncture  points. 

After  all  these  functions  have  been  completed,  GRIDIN  turns  over  control  of  the 
program  to  subroutine  SGRID2.  After  cyclmg  through  subroutines  SGRID2  to  SGRID8, 
control  returns  to  GRIDIN,  which  then  stores  the  plot  under  its  identification  number 
in  the  data  set  BLADE.  PLOTS,  by  calling  subroutine  ENDID.  The  plot  is  thus  saved 
for  redisplay  at  some  future  time.  GRIDIN  does  display  the  plot  at  this  time  as  well, 
on  the  device  previously  selected  by  the  user.  The  slice  counter  ISLICE  is  then  incre¬ 
mented  by  1,  and  GRIDIN  loops  back  to  read  the  TACTGRID  input  for  the  next  blade 
slice.  The  TACTGRID  program  is  terminated  by  GRIDIN,  upon  exhaustion  of  the  input 
data,  with  a  conditional  READ  statement. 


Subroutine  TIME 

Subroutine  TIME  is  a  general-utility  subroutine  available  on  the  IBM  TSS/360-67 
computer  at  the  Lewis  Research  Center.  It  places  the  numei’ical  characters  for  the 
month,  day,  year,  hour,  and  minute,  when  called,  into  array  ALPH.  This  string  of 
characters  is  printed  on  the  first  page  of  the  output  and  displayed  in  the  upper  right 
corner  of  the  plot  (fig.  4). 


Subroutine  TFORM 

TFORM  is  the  main  program  of  the  code  in  reference  6 .  This  main  program  is 
converted  into  a  subroutine  for  use  by  TACTGRID.  As  in  reference  6,  TFORM  uses 
subroutines  ROOT,  FUNCT,  and  BLCD;  but  the  spline  curve  subroutine  of  reference  6, 
SPLN22,  has  been  replaced  by  an  improved  spline  program,  SPLISL.  The  function  of 
TFORM  and  its  slave  subroutines  is  described  in  detail  in  reference  6. 


Subroutine  SPLEL 

Subroutine  SPLEL  calculates  the  first  and  second  derivatives  of  a  cubic  spline 
curve  at  the  spline  points.  It  is  found  in  the  MERIDL  code  and  described  in  refer¬ 
ence  8.  The  SPLEL  routine  is  based  on  the  end  condition  of  specified  slopes  at  the 
two  end  points. 
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Subroutine  BEG  ID 


Subroutine  BEGID  is  a  Lewis  Research  Center  graphics  system  subroutine.  It 
marks  the  beginning  of  a  set  of  common  plot  orders  and  assigns  an  identification  num¬ 
ber  to  the  set.  By  referring  to  the  Lewis  graphics  system  manual  (ref.  9),  the  user 
should  be  able  to  replace  the  graphics  subroutines  in  TACTGRID  by  graphics  subrou¬ 
tines  of  comparable  function  available  at  his  facility  The  word  descriptions  of  these 
subroutines  are  taken  from  reference  9. 


Subroutine  SCALE 

Subroutine  SCALE  is  a  Lewis  graphics  system  subroutine.  It  scans  an  array  of 
plot  coordinates  to  determine  maximum  and  minimum  values.  These  values  are  used 
by  the  graphics  package  to  set  the  scale  of  the  plot. 


Subroutine  SETFRM 

Subroutine  SETFRM  is  a  Lewis  graphics  system  subroutine.  It  allows  the  user  to 
specify  the  number  of  frames  in  his  plot,  creating  a  larger  physical  plotting  area  in  the 
x-direction. 


Subroutine  NAXIS 

Subroutine  NAXB  is  a  Lewis  graphics  system  subroutine.  It  enables  the  user  to 
specify  display  boundaries  when  no  axes  are  to  be  drawn.  GRIDIN  uses  this  subroutine 
to  position  the  blade  slice  plot  within  the  frame . 


Subroutine  SCISS 

Subroutine  SCISS  is  a  Lewis  graphics  system  subroutine.  It  allows  the  user  to  set 
image  scissoring  options.  Image  scissoring  involves  truncating  those  portions  of  a 
display  unit  that  extend  beyond  the  boundaries  of  the  plot  or  device.  GRIDIN  uses  this 
subroutine  to  set  the  scissoring  at  screen  boundaries  but  employs  an  option  which  al¬ 
lows  continued  point  examination  after  scissoring. 
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Subroutine  GPLOT 


Subroutine  GPLOT  is  a  Lewis  graphics  system  subroutine.  It  enables  the  user  to 
plot  single  or  multiple  curves  or  lines  as  vector,  point,  symbol,  or  vector-symbol 
plots , 


Subroutine  ENDID 

Subroutine  ENDID  is  a  Lewis  graphics  system  subroutine.  It  indicates  the  end  of 
a  set  of  common  plot  orders  assigned  an  identification  number  by  subroutine  BEGID. 

At  the  user' s  option,  this  subroutine  will  save  the  plot  order  set  in  a  user-provided 
partitioned  data  set  (VPAM)  or  in  virtual  memory.  GRIDIN  uses  subroutine  ENDID  to 
save  the  plots  of  the  individual  blade  slices  in  the  VPAM  data  set  BLADE .  PLOTS.  The 
identification  number  of  each  slice  plot  is  the  same  as  the  slice  number  ISLICE. 


Subroutine  DISPLA 

Subroutine  DISPLA  is  a  Lewis  graphics  system  subroutine.  It  defines  the  end  of  a 
display  or  plot  and  initiates  transmission  of  plot  orders  to  the  plotting  device. 


Subroutine  TERM 

Subroutine  TERM  is  a  Lewis  graphics  system  subroutine.  It  closes  the  graphics 
data  set  and  performs  cleanup  functions.  GRIDIN  uses  subroutine  TERM  to  terminate 
the  plottii^  and  close  the  data  set  BLADE .  PLOTS  before  terminating  the  program  it¬ 
self. 


Subroutine  SGRID2 

Subroutine  SGR1D2  determines  those  TFORM -calculated  points  that  are  in  the 
spline  portion  of  the  blade.  It  calculates  the  surface  lengths  of  the  suction-  and  pres¬ 
sure-side  outer  walls  (STPNEW  and  SBTNEW)  in  order  to  determine  where  to  begin 
tapering  the  wall  and  channel  thicknesses  in  accordance  with  the  input  arrays  THKIT, 
THK2T,  THK3T,  THKIB,  THK2B,  and  THK3B  described  previously.  SGRID2  plots  the 
outside-wall  spline  contours  for  debugging  purposes.  Subroutine  SGRID2  plots  the  ra¬ 
dial  lines  defining  the  ends  of  the  circular  regions  on  the  suction  and  pressure  sides, 
for  reference  on  the  blade  slice  plot.  For  slices  beyond  the  first,  these  lines  need  not 
coincide  with  the  last  circular- region  stations. 
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SGRID2  calculates  all  the  nodal  points  (and  intermediate  vector  plotting  points)  in 
the  circular  region  of  the  blade  shell.  It  plots  the  circular- region  horizontal  and  ver¬ 
tical  orthogonals  in  the  blade  shell. 

Subroutine  SGRID2  calculates  the  normal-direction  and  tangential-direction  arc 
length  distances  (arrays  NCB,  NCT,  SCB,  and  SCT  of  the  output)  along  the  grid  or¬ 
thogonals  in  the  blade  shell  for  the  circular  portion  of  the  blade. 


Subroutine  SPINSL 

Subroutine  SPINSL  is  a  cubic  spline-curve  program  used  for  interpolation,  includ¬ 
ing  interpolation  of  first  and  second  derivatives.  It  is  found  in  the  MERIDL  code  and 
described  in  reference  8.  It  uses  end  conditions  of  specified  end-point  slope.  SPENSL 
is  an  alternate  entry  point. 


Subroutine  SGRID3 

Subroutine  SGRID3  calculates  the  ariays  of  points  on  the  horizontal  orthogonals  for 
the  suction  and  pressure  sides  in  the  blade  walls.  The  first  horizontal  orthogonal  is 
the  blade  outer  surface .  The  fourth  horizontal  orthogonal  is  the  blade  iimer  surface . 
SGRID3  uses  the  root-finding  subroutine  ROOT  in  two  levels  of  iteration  (see  section 
PROGRAM  PROCEDURE  FOR  FIRST  SLICE)  to  cut  the  straight  trailing-edge  channel 
and  employs  subroutine  INRSCT  to  locate  the  intersections  of  the  outer  surface  normals 
with  the  channel  walls .  These  points  are  used  to  fit  the  spline  curve  in  region  VI  of  the 
blade  (fig.  10). 

Subroutine  SGRID3  constructs  the  three  horizontal  orthogonals  for  the  impingement 
crossflow  channel  and  insert  walls.  The  second  horizontal  defines  the  insert  outer  con¬ 
tour.  The  third  horizontal  orthogonal  defines  the  insert  inner  contour. 

Subroutine  SGRID3  calculates  the  surface  lengths  of  the  blade-shell  iimer  contours 
for  use  in  determining  the  distribution  of  station  positions,  as  described  previously. 
SGRID3  plots  the  blade-shell  inner  surface  contours  for  debugging  purposes. 

Subroutine  ROOTl 

Subroutine  ROOTl  is  the  bisection-method,  root-finding  subroutine  ROOT  in  the 
older  version  of  the  MERIDL  code  (ref.  10).  Two  copies  of  subroutine  ROOT,  one 
called  ROOTl  and  the  other  ROOT2,  are  used  by  subroutine  SGRID3.  The  two  names 
are  necessary  because  two  nested  loops  use  the  bisection  procedure  for  iteratively 
determining  the  orientation  of  the  trailing-edge  channel  and  a  subroutine  may  not  call 
itself  (fig.  21). 
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Subroutine  FUNCl 


Subroutine  FUNCl  is  the  external  subroutine  in  the  argument  of  ROOTl  that  calcu¬ 
lates  the  function  whose  root  is  to  be  found.  In  this  case,  for  the  outer  iteration  loop 
(see  section  PROGRAM  PROCEDURE  FOR  FIRST  SLICE),  this  function  is  the  directed 
distance  sgn(Xj  -  Xj)*|lJ| /THKHOL  of  figure  16(a),  where  sgn(a)  =  1  for  a  20  and 
sgn(a)  =  -1  for  a  <  0. 


Subroutine  ROOT2 

Subroutine  RDOT2  is  the  second  copy  of  subroutine  ROOT,  referred  to  previously, 
which  is  called  by  subroutine  FUNCl  for  the  inner- loop  iteration  (fig.  21). 

Subroutine  FUNC2 

Subroutine  FUNC2  is  the  external  subroutine  in  the  argument  of  ROOT2  that  calcu¬ 
lates  the  function  whose  root  is  to  be  found.  In  this  case,  for  the  inner  iteration  loop 
(see  section  PROGRAM  PROCEDURE  FOR  FIRST  SLICE),  this  function  is  the  differ¬ 
ence  in  lengths  (lEB  I  -  IFD  I  )/THKHOL  of  figure  15(a). 

Subroutine  INWSCT 

Subroutine  INWSCT  is  a  program  for  finding  the  intersection  of  two  spline  curves 
by  the  method  of  slopes.  It  is  the  INRSCT  subroutine  of  reference  8  with  two  changes. 
The  first  change  is  that  the  search  begins  at  the  right  end  point  of  the  search  interval 
rather  than  in  the  middle .  The  other  change  is  that  the  first  curve  is  fit  with  subrou¬ 
tine  SPINSL  and  the  second  curve  with  subroutine  SPLINT,  rather  than  both  curves  be¬ 
ing  fit  with  SPLINT,  as  is  done  in  INRSCT.  Subroutine  INRSCT  is  described  in  refer¬ 
ence  8 .  INWSCT  is  used  by  FUNC2  for  finding  the  intersection  of  the  trailing-edge- 
channel  bottom  surface  with  the  blade  pressure-side  inner  surface. 


Subroutine  SPLINT 

Subroutine  SPLINT  is  a  cubic  spline-curve  program  used  for  interpolation,  includ¬ 
ing  interpolation  of  first  and  second  derivatives.  It  is  found  in  the  MERIDL  code  and 
described  in  reference  8 .  It  uses  the  end-point  condition  that  the  second  derivative  at 
either  end  point  is  one-half  that  of  the  next  spline  point.  SPLENT  is  an  alternate  entry 
point. 
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Subroutine  INRSCT 


Subroutine  INBSCT  is  the  intersection-finding  program  (mentioned  previously)  of 
reference  8.  It  is  described  in  that  reference. 


Subroutine  SPLNSL 

Subroutine  SPLNSL  is  a  cubic  spline-curve  program  used  for  interpolation,  includ¬ 
ing  interpolation  of  first  and  second  derivatives.  It  is  a  hybrid  of  the  two  spline- curve 
programs  SPLINT  and  SPINSL  of  reference  8.  SPINSL  is  a  spline-curve  program  that 
requires  specified  slopes  for  its  end  condition.  SPLINT  does  not  require  the  specifica¬ 
tion  of  slopes  for  its  end  condition.  SPLNSL  uses  a  specified  slope  for  its  left  end  con¬ 
dition  and  the  SPLINT  end  condition  for  its  right  end  condition.  Subroutine  SPLNSL  is 
used  for  spline  curve  fitting  of  the  midwall  and  impingement  crossflow-channel  horizon¬ 
tal  contours,  where  the  left-end-point  slope  is  known  (tangential  to  a  circular- region 
contour)  but  the  right-end-point  slope  is  not.  SPONSL  is  an  alternate  entry  point. 


Subroutine  SGRID4 

Subroutine  SGRID4  calculates  the  horizontal  orthogonals  (approximate  streamlines) 
for  the  trailing-edge  flow  channel  (regions  V  and  VI  of  fig.  10).  It  locates  the  position 
of  the  plenum  end  within  the  blade.  The  subroutine  determines  station  positions  in  the 
blade  shell  and  generates  the  vertical  orthogonals  for  the  blade  shell  and  in  the  trailing- 
edfee  flow  channel. 

Subroutine  SGRID4  calculates  the  normal-direction  arc  length  distances  (arrays 
NBB  and  NTT  of  the  output)  along  the  vertical  orthogonals  in  the  blade  shell  for  the 
spline  portion  of  the  blade.  If  needed,  SGRID4  redefines  the  midwall  horizontal  orthog¬ 
onal  and  reconstructs  the  vertical  orthogonals  in  the  blade  shell,  iterating  until  the  ver¬ 
tical  orthogonals  are  split  in  two . 

Subroutine  SGRID4  plots  the  blade-shell  vertical  orthogonals  in  the  spline  region  of 
the  blade. 

The  procedures  employed  in  subroutine  SGR1D4  are  discussed  in  the  section  PRO¬ 
GRAM  PROCEDURE  FOR  FIRST  SLICE. 


Subroutine  INTSCT 


Subroutine  INTSCT  is  identical  to  subroutine  INRSCT  except  that  the  first  curve 
is  fit  with  subroutine  SPINSL  rather  than  subroutine  SPLINT.  It  is  used  instead  of 
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INRSCT  in  the  orthogonalization  procedure  when  the  horizontal  orthogonal  has  known 
end-point  slopes. 


Subroutine  FUNC4 

Subroutine  FUNC4  is  the  external  subroutine  in  the  argument  of  a  RCXDTl  call  in 
subroutine  SGRID4.  It  is  used  in  determining  the  intersection  of  the  suction-side  ple¬ 
num  external  contour  with  the  pressure-side  plenum  external  contour.  This  intersec¬ 
tion  point  then  establishes  the  search  interval  in  which  to  seek  the  location  of  the  ple¬ 
num  end.  FUNC4  calculates  the  difference  |yj^(x)  -  Y2(x)] /sDUMP,  where  is  the 
ordinate  of  a  point  on  the  suction-side  contour  and  Y2  is  the  corresponding  ordinate 
on  the  pressure-side  contour.  Of  course,  when  Yj^(x)  =  Yglx),  the  intersection  has 
been  found. 


Subroutine  FUNC3 

Subroutine  FUNC3  is  the  external  subroutine  in  the  argument  of  a  ROOTl  call  in 
subroutine  SGRUM.  It  constructs  a  vertical  orthogonal  link  (by  the  method  of  ref.  5) 
between  the  pressure-  and  suction-side  plenum  external  contours,  from  a  search  point 
on  the  pressure-side  contour.  It  calculates  the  length  of  this  vertical  orthogonal  link, 
ENDLNG .  FUNC3  then  evaluates  the  function  (SDUMP  -  ENDLNG)/SDUMP,  whose  root 
is  found  by  BCXDTl,  to  determine  the  position  of  the  plenum  end. 


Subroutine  INMSCT 

Subroutine  INMSCT  is  the  same  as  subroutine  INTSCT  except  that  the  first  curve 
is  fit  with  subroutine  SPLNSL  rather  than  with  subroutine  SPINSL.  It  is  used  instead 
of  INBSCT  in  the  orthogonalization  procedure  when  the  horizontal  orthogonal  has  a 
known  left  end-point  slope. 


Subroutine  SGR1D5 

Subroutine  SGRID5  calculates  the  tangential-direction  arc  length  distances  (arrays 
SBB  and  STT  of  the  output)  along  the  horizontal  orthogonals  in  the  blade  shell  for  the 
spline  region  of  the  blade.  It  interpolates  points  between  grid  points  in  the  spline  re¬ 
gion  of  the  blade  shell  in  order  to  plot  the  horizontal  orthogonals  and  to  calculate  arc 
lengths  between  grid  points  accurately .  The  subroutine  plots  the  horizontal  orthogo¬ 
nals  in  the  blade  shell  for  the  spline  region  of  the  blade . 
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Subroutine  SGRID6 


Subroutine  SGRID6  constructs  the  orthogonal  grid  for  the  impingement  crossflow 
channel  by  generating  vertical  orthogonals  from  the  blade  inner-surface  nodal  points. 
The  routine  calculates  the  normal-direction  arc  length  distances  (arrays  NBB  and  NTT 
of  the  output)  along  vertical  orthogonals  in  the  impingement  crossflow  channel. 

If  needed,  SGRID6  redefines  the  mean  streamline  and  reconstructs  the  vertical  or¬ 
thogonals  in  the  impingement  crossflow  channel  -  iterating  until  the  vertical  orthogo¬ 
nals  are  split  in  two.  SGRID6  plots  the  impingement-crossflow-channel  vertical  or¬ 
thogonals  in  the  spline  region  of  the  blade . 


Subroutine  SGRID7 

Subroutine  SGRID7  calculates  the  tangential-direction  arc  length  distances  (arrays 
SBB  and  STT  of  the  output)  along  the  horizontal  orthogonals  in  the  impingement  cross- 
flow  channel  for  the  spline  region  of  the  blade.  It  interpolates  points  between  grid 
points  (and  intersection  points  of  vertical  orthogonals  with  the  plenum  inner  and  outer 
contours)  in  the  spline  region  of  the  impingement  crossflow  channel  in  order  to  plot  the 
horizontal  contours  and  to  calculate  arc  lengths  between  stations  accurately .  This  sub¬ 
routine  plots  the  horizontal  contours  in  the  channel  (mean  streamline  and  insert  inner 
and  outer  surfaces)  for  the  spline  region  of  the  blade. 

Subroutine  SGRID7  calculates  the  plenum  perimeter,  cross-sectional  area,  and 
hydraulic  diameter  and  plots  the  insert  end  length.  This  is  the  dump  hole  opening  if  a 
"cutoff"  design  for  the  insert  is  used.  SGRID7  calculates  all  the  nodal  points  (and  in¬ 
termediate  vector  plotting  points)  in  the  circular  region  of  the  impingement  crossflow 
channel.  It  plots  the  circular-region  horizontal  and  vertical  orthogonals  in  the  channel. 
SGRID7  calculates  the  normal-direction  and  tangential-direction  arc  length  distances 
(arrays  NCB,  NCT,  SCB,  and  SCT  of  the  output)  along  the  grid  orthogonals  in  the  im¬ 
pingement  crossflow  channel  (including  insert  surfaces)  for  the  circular  portion  of  the 
blade.  SGRID7  calculates  the  channel  he^hts  and  channel  lengths  for  the  stations  in 
the  traUing-edge  chaimel,  beginning  with  the  mixing  station.  It  plots  ibe  mean  stream¬ 
line  and  channel  vertical  orthogonals.  SGRID7  prepares  the  TACTl  input  data  set  for 
the  individual  blade  slice  (fig.  13(f)). 


Subroutine  SGRID8 

Subroutine  SGRID8  calculates  the  arrays  of  points  used  in  constructing  the  arrows 
to  the  calculation  stations  and  plots  these  arrows  with  the  station  numbers  indicated. 
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SGRID8  calculates  the  scale  of  the  plot  and  prints  the  scale  and  legend  on  the  plot,  as 
well  as  all  other  plot  labels.  SGRID8  transforms  the  output  coordinates  of  the  individ¬ 
ual  blade  slices  to  the  TFORM  coordinate  system  set  up  for  the  first  slice,  which  is 
taken  as  global.  The  subroutine  then  prints  out  these  transformed  coordinates  in  the 
form  depicted  in  figure  13(g). 


Subroutine  NUMBER 

Subroutine  NUMBER  is  a  Lewis  graphics  system  subroutine.  It  allows  the  user  to 
convert  an  integer  or  real  number  to  an  array  of  printable  characters.  It  is  used  by 
subroutine  SGRID8  to  convert  the  station  numbers  to  character  data  (A4  format)  so  that 
they  can  be  printed  on  the  blade  sUce  plot  by  subroutine  CHARS . 


Subroutine  CHARS 

Subroutine  CHARS  is  a  Lewis  graphics  system  subroutine.  It  allows  the  user  to 
print  character  data  anywhere  on  a  plot.  It  is  used  by  subroutine  SGRID8  for  printing 
the  labels  and  station  numbers  on  the  plots  of  the  blade  slices . 
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Figure  1.  -  Three-dimensional  schematic  view  of  a  turbine  blade  equipped  with  impingement 
insert  and  pin  fins. 


Impingement  ^  ^  Coolant  channel  (fed 

/  '  by  impingement) 


Detail  of  nodes 

Figure  2.  -  Blade  geometric  model. 


,4  16  18 


0  .1 

in. 

Figure  3.  -  Cross  section  of  sample  problem  blade. 


«  .t  l*(C“ 

sc*; f  :s  6. :'s»  ' !■{ s  s; /f 


CP055  SECTION  OF  SLICE  Nu»*BEP  ! 

(a)  Hub  section  of  sample  blade  problem. 
Figure  4.  -  Scaled  computer  plots  of  blade  slices. 
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BOPOSED  first  stage  BQTQR,  MEAN 


THE  DATE  OF  THIS  BUN  IS  l?/06/^0  18:05 


«  I  2  9  NN 


0  .1  INCH 

sent  IS  ♦.2J/3  TIMC5  5l2t 


CBOSS  SECTION  OF  SLICE  NUMBEB  2 

(b)  Mean  section  of  sample  blade  problem. 
Figured  -  Continued. 
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I  THSPl(-) 


Figure  7.  -  Transformed  coordinates  and  transformation  constants. 
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TACTGRID 

input 

y.o  7.0 

3  2  1 

.  333  3  33  333.  ec 666666 6. 333333  333.66t>(;(i66 66 

4  6  S  S 

6  5  5  10  1 

6GE0  THK  1T*51*.010,THK2T=51*.03S,THK31=S1*.03b,ldKlB  =  !)1*. 
.0001  ,C001  .030  .010  .080  8.875 

0.0  0.0 

0  10,THK2E  =  5U.  03  5,THK3b-S1  •.  035,OEbi; 
8.50 

PBOPOSED 

FIRST  STAGE  ROTOR,  HUE 

0.03429 

-0,06000  0.21589 

39.37009 

0. C5000 

0. 003046 

0.0006699  52.00 

-56.28 

7. CO 

IrUK/Vi 

0.00000 

0.003810  0.01143 

0. 01788 

0-02486  0.02921 

O.OOOOC 

input 

0.00000 

0.02082  0.03001 

0.C23C7 

O.COIhoO  -0.01812 

!  O.OOOOC 

0.003048 

0.0008899  22.00 

-40.28 

5.00 

O.OOOOC 

0.01016  0.02032 

0.03046 

C. COCOC 

0.00000 

-0.006999  -0.01S0C 

-0.04393 

O.OOOCO 

TACT- 

63E0  THKlT=51*.OlO,THK2r=51* 

.035,THK3T 

=514.035,THK1B=51*. 

0  1&,lHK2B-51*.C35,:fiK3c  =  5  1*.  035,6FND 

.0001 

.0001  .030 

.010 

.070  9.625 

8.875 

GRID 

0.0 

0.0 

priOPOSEc 

FIRST  STAGE  RCTOP,  .1EAN 

0.03429 

-0.05405  0.2349499 

39-37009 

0.C5C00 

0,003048 

0.0006899  42.00 

-58.58 

7.00 

TFORM 

0.00000 

0.003610  0.C1143 

0.01786 

0.02409  C.C7921 

C.C^CCC 

0.00000 

0.C17S1  0.C2487 

0.01850 

0.0009950  -0.02175 

O.OJCOt 

0.003046 

0.0006899  12.00 

-42.58 

5.CC 

0. 00000 

0.01016  0.02032 

0. 03048 

O.COOOO 

0.00000 

-0.009140  -0.01738 

-0.04777 

O.OOCOO 

TACT- 

60EC  THK1T*51*.010,THK2T*51* 

.035,THK3T 

a51*.03S,THK1B*51*. 

010,'!HK*.E  =  51*.  055,THK.3e*51*.  J35,&ESD 

.0001 

.0001  .030 

.010 

.060  1C. CO 

9.625 

GRID 

0.0 

0.0 

PROPOSED 

FIRST  STAGE  RCTOt,  T 

IP 

0.03429 

-0.05750  0.254C 

.39.37009 

0. C5000 

C.C03048 

0,0006899  32.00 

-60. tb 

7.00 

TFORM 

0.00000 

0.003810  0.01143 

0.01788 

0. 02340  0.02921 

o.coccc 

O.COOOO 

0.01470  0.02051 

0.01461 

0.00O8370  -0.02444 

I^-OOCOO 

input 

0.003048 

0.0008899  2.0C 

-44.68 

5.  CO 

O.OGOOO 

0.01016  0.02C3Z 

0.03046 

O.COOOO 

0.00000 

-0.01095  -0-C195C 

-0.  05105 

0.00000 

(a)  Sample  problem  input  data  set  listing. 
Figure  12.  -  TACTGR ID  input  format. 
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THE  DATE  OP 

THIS  RUN  IS  08/29/79 

VRO 

HBU 

9.000 

7.000 

NSLICE 

OU  NPEAHE 

3 

2  1 

PIT 

F2T  FIB 

C.333 

0.667  0.333 

NPT 

NIT  N2T 

4 

6  5 

NIB 

N2B  N36 

6 

5  5 

6CKANLS  »SLICE=  3  ,aST»=69,  I'lKItSi 

09:37 


F2B 

0.667 

H3T 

S 

NUB 

10 

2«£E!iD 


N5B 

5 


(a)  TACTGRID  and  TACTl  overall  input. 


I 

IHK1T  (I) 

TtiK2T  (I) 

1 

0.0100 

0.0350 

2 

0.O1CO 

0.0350 

3 

0.0100 

0.0350 

4 

0.0100 

0.0350 

5 

0.0100 

0.0350 

6 

0.0100 

0.0350 

7 

0.0100 

0.0350 

e 

0.0100 

0.0330 

9 

0.0100 

0.0350 

10 

0.0100 

0.0350 

11 

0.0100 

0.0350 

12 

0.0100 

0.0350 

13 

0.0100 

0.0350 

14 

0.01G0 

0.0350 

15 

0.0100 

0.0350 

16 

0.0100 

0.0350 

17 

0.0100 

0.0350 

1b 

0.0100 

0.0350 

19 

0.0100 

0.0350 

20 

0.0100 

0.0350 

21 

0.0100 

0.0350 

22 

0.0100 

0.0350 

23 

0.0100 

0.0350 

24 

0.0100 

0.0350 

25 

0.0100 

0.0350 

26 

0.0100 

0.0350 

27 

0.0100 

0.0350 

28 

0.0100 

0.0350 

29 

0.0100 

0.0350 

3C 

0.C100 

0.0350 

31 

0. 0100 

0.0350 

32 

0.0100 

0.0350 

33 

0.0100 

0.0350 

34 

0.0100 

0.0350 

35 

0.0100 

0.0350 

36 

0.0100 

0.0350 

37 

0.0100 

0.0350 

36 

0.0100 

0.0350 

39 

0.0100 

0.0350 

40 

0.0100 

0.0350 

41 

0.0100 

0.0350 

42 

0.0100 

0.03S0 

43 

0.0100 

0.0350 

44 

0.0100 

0.0350 

45 

0-0100 

0.0350 

46 

0.0100 

0.0350 

47 

0.0100 

0,0350 

46 

0.0100 

0.0350 

49 

0.0100 

0.0350 

50 

0.0100 

0.0350 

51 

0-0100 

0.0350 

THK(1) 

THK  (3) 

THKHOL 

0.000 

0.000 

0.030 

Z 

THETA 

0.000 

0.000 

THK3T(I) 

ThKlP(l) 

THK2B(I) 

0.0350 

C.0100 

0.0350 

0.0350 

0.O1CO 

0.0350 

0.0350 

0.0100 

0.0350 

0.035C 

0.010C 

0.0350 

0.0350 

G.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

:.oioo 

0. 0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

C.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

c.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0. 0350 

0.0350 

0.0100 

0.0350 

0.0350 

c.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

c.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.C10C 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

c.0100 

0.0350 

0.C350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

c.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

C-0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0-0350 

0.0350 

0.0100 

0.0350 

0.0350 

c.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

c.0100 

0.0350 

C.0350 

0.0100 

0.0350 

0.0350 

0-0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

C-0100 

0. 0350 

0.0350 

0-01C0 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

0-0350 

0-0100 

0.0350 

0.0350 

0.0100 

0.0350 

0.0350 

0.0100 

0.0350 

THKHAL 

5D0HP 

BBO 

0.010 

0.080 

8.875 

(bl  TACTGRID  geometrical  input  for  blade  slice. 
Figure  13.  -  TACTGRID  sample  output. 


THK3c  (1) 
0.0350 
0.0550 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.035C 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.C350 
0.0350 
0.03SO 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
C.0350 
0.0350 
PHI 
6.500 
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PSOPOSED  PIBST  STAGE  BOTOB,  HUB 


CHORD 

S16R 

RHI 

SCALE 

DELX 

0.3429000E-01 

-0.5000000E-01 

0.2156900 

39,37009 

O.SOOOOOOE-OI 

BLADE  SURFACE 

1  —  OFFER  SURFACE 

B11 

R01 

BETH 

BET01 

SPLROI 

0.  J04eO0OE-02 

0.8898999E-03 

52.00000 

-56.28000 

7.O0C0O0 

nSPI  ABRAY 

0.0000000 

0.3810000E-02 

0.  11430 OOe-OY 

0, 1788000E-01 

O.248t)OOOE-01 

0.2921000E-01 

0.0000000 

THSP1  ABRAY 

0.0000000 

0.2082000E-01 

0. 3001000E-01 

0. 2307000 E-01 

0- 1460000E-02 

-0.1812000E-01 

O.OOOCOOO 

BLADE  SURFACE 

2  L04EB  SURFACE 

EI2 

B02 

BETI2 

BET02 

SPLUOP 

0.3048000E-02 

0.8898999B'03 

22.00000 

-40.28000 

5.000000 

NSP2  ARRAY 

O.OQQQOOO 

0.1016QQ0E-01 

0.2032QQQE-01 

Q.304800QE-01 

O.OOOOOOQ 

THSP2  ARRAY 

0.0000000 

-0.6998997E-02 

-0.  1500000E-01 

-0.43930008-01 

0.0000000 

(c)  TFORM  input  for  blade  slice. 


BLADE  DATA  AT  INPUT  SPLINE  POINTS 
BLADE  SUBFACE  1 


n 

THETA 

DERIVATIVE 

2N0  DERIV. 

0.64614E-03 

0.86921E-02 

5.9287 

-1810.4 

0. 38100E-02 

0.20820E-01 

2.5064 

-352.98 

0.  11430E-01 

0. 30010E-01 

-0.49750E-01 

-317.91 

0.  17880E-01 

0.23070E-01 

-2.1031 

-318.79 

0.24B60E-01 

0.  14800E-02 

-3.9605 

-213.40 

0-29210B-01 

-0.  iai20E-01 

-5.1321 

-325.26 

0.34140E-C1 

-0.47712E-01 

-6.9401 

-408. 18 

ELADE  SURFACE  2 

H 

THETA 

DERIVATIVE 

2ND  DERIV. 

0.4ld98E-02 

-0.13090E-C1 

1.8714 

-321.57 

0-  10160E-01 

-0.69990E-02 

0.27786 

-212.27 

0.20320E-01 

-0.  lSOOOE-01 

-1.8399 

-204.60 

0. 304eOE-01 

-0.a3930E-01 

-3.8232 

-185.82 

0. 32825E-01 

-0.53145E-01 

-3.9254 

96.618 

(d)  TFORM  profile  smoothness  check. 


NO.  OF  POINTS  * 

30  PHI  s 

-15.74  CIGRBES 

X 

Y  LONER 

Y  UPPER 

0.00000 

0.12000 

0.12000 

0.50000E-01 

0.22533E-01 

0.24Q55 

0.10000E  00 

0.167d6E-02 

0.29761 

0.  15000 

0.381C2E-02 

0.34204 

0.20000 

0.30113E-Q1 

0.37663 

0.250C0 

0.64618E-01 

0.4062S 

0.30000 

0.93299E-01 

0.43177 

0.35000 

0. 11707 

0.44969 

0.40000 

0. 13670 

0.46259 

0.45000 

0. 15277 

0.47C86 

0.50000 

0.16556 

0.47482 

0.55000 

0.17529 

0.47476 

0.60000 

0.18207 

0.47094 

0.65000 

0. 18607 

0,46356 

0.70000 

0.18740 

0,  45297 

0.75000 

0.16622 

0.43944 

0.80000 

0.18261 

0.42331 

0.85000 

0.17671 

0.4C485 

0.900GO 

0.16860 

0. 36428 

0-95000 

0.15840 

0,36182 

1.00000 

0.14621 

0,33759 

1.0500 

0.13211 

0.31146 

1.1000 

0.11619 

0.26335 

1. 1500 

Q.98519E-01 

0.25310 

1.2000 

0.79161E-01 

0.22065 

1.2500 

0.58239E-01 

0.18596 

1.3000 

0.35e4lE'01 

0.14909 

1.3500 

0.127621-01 

0.10997 

1.4000 

0.31670E-02 

O,68650E-O1 

1.4205 

0.35035E-01 

0.35035E-01 

(e)  TFORM  output. 


Figure  13.  -  Continued. 
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sShi)  :2TA  =  U,  I;;TB  =  1  2,THK--0.  CIOGO,  0.  0  3  S  0  0,  0  .  U  3  SCO,  1  DLX  =  0 . 052d  9 , 0.  05  1  76 , 0  .  C  4  9  76 , 0  .  G477c3,G.04S/a,tfNr. 
••  ;eg  I  STA  =  1  3, 131b=  1  3, THK  =  0. 01  CGG,  0-0  3S00,  0, 03500,1  CLX=0.03337,0.  03364,0.  03466,0.  03546,  G.iJJt  29,  lend 
c-.;r:(i  IS!  A  =  1  4,  lSTD  =  l4,l!iK-^0.0U00,0.O^S0O,0. 0  3  500,  *1  DLX  =  0.  05090, 0.  050  2  1 ,  C  .  0  4  9  OG  ,  0  .  G4  77  c  ,  G  -  o4  c  S  7  ,  l  h  N  [, 
GG.mJ  ISTA  -  is,  IGTD  =  1S,THK=0.()tL00,  C-0  3..UG  ,0.0  3S00,1CLX  =  0.0314  1  ,u.  jJ  jb7,a.  03406,0.  J  jL*4o,  G,  v3f.zs,  6  t : 

GF.C  ISTA=--1  6, 1016=^16  ,THK--O.O10C0, 1-.  03SJ0,0. 03500,1  CLX  =  C.051  2  4,C-05c47,C.C49l3,0.C4  776,G.:4^44,» 
6GFL'  iSTA=17, 1ST E=  17, THK-O. 01000,  0.03500,0.  Ci^'OC, 7  LLX  =  0. 0  3 34  7,0.  0 3  392 ,  C .  o 3 4 70 , 0 .  G  354 c  ,  G-  G 3 c 2 c  ,  c  x  . 
{■i;EO  I0TA^16,IETB  =  16,THK  =  0.01CCG,0.0350O,U.Uj‘'CO,7CLX  =  0.051S0,O.  D50b7,0.04y2j,O.U477o,0.04bi3,c  i:Sj 
eCiU  r.;iA  -  19, 1676=  19,  rHK=C.  CUOO,  G.0>5CG,0.0  3*^00,  TDLX=0,  0335  s,  U.  G33yo,C.Ci473,C.03548,C.. 
f-Gt'i.)  .STA  =  -0,  ISTB=20,THK  =  0.  0  1000,  0.  0  3  500,  0 , 0  3  SCO,  7  DL  X=0.  0C2  On  ,  o  .  o  6  1  ^  3 , 0 .  ;5  5l6,G.GS7j4,G.i-Jj4  9,t£;i2 
t-GEO  ISTA*  2  1,ISTB=2  1  ,THK  =  0.  01000,0.03500,0. 0?SGG,1  CLX  =  C.  040  39, 0.  L4C6fe,  0.  04  1  72,0. 0425  7,  G..  G4  J4.,i  i..j 
GGEO  I5TA=  22,  ISTB  =  2’2,  THK=0.  01000,  0.03  500,0.03500,1  DLX  =  0.C£  22  5,0.  06  1  16,  C.  059  i.s,C.  JS7  34, 0..  ^  o'  4  5, 6  £  NL 
SGSO  ISTA  =  2  3,ISTB  =  2  3,1HK  =  0.01GGU,  0.03  50C  ,0.  U3500,TCLX-0.  0404  6,0.  04C93, O.G4l7S,u.  04  257, j.:,4:4v,6E.‘.j 
6GEO  ISTA=24,  ISTB  =24,  TBK  =  0.01  000, 0.03  500,0. 03  5CC,lDLy =0.0622  3,0. 06  1l4,C.C5  9  24,C.C5734,C.i5S-4  3,^ESr 
5GEO  ISTA=25,ISTB  =  25,THK  =  0.01000,  0. 03500, 0 . 0  350C,7DLX  =  0. 0404  7 , 0-  C4  0  94,0.04l7D,0.04257,i^-G4>j‘,t£NL 
SGEO  IGTA  =  2  6, ISTB  =  26,THK  =  U.01 000, 0, 03600,0 . 03500, 7  CLX=0.0(2 1 2, 0.  06 1 06 , C. G59 20 , 0, OS? 34 , C. ^ S Si c ,  i  £  NO 
&GlO  I‘iTA  =  2  7,  :STB=2  7,THK  =  0.01000, 0.  03  500,0,Oj50C,  TDLX=O.OU04  8,C.04C95,O.C-«17G,O.G4  2  57,J.C433y,6£'iL 
5  GEO  I5TA  =  28,ISTB=28,7HK  =  0.OlC0O,G.O35OO,O-OJsOO,  7CLX=0.061B?,'',  OtCco  ,0.  059  1'',C.  0S7Ja,  0.  0  S  357,  G  i.:;2 
flGEr>  ISrA  =  2y,  IsrB=2  9,  THK=G.  01OGO,  0.  03500,  0. 0  3 ‘•oG,  TDIX=0.  0405  1,  C.  04  C57,  0.  v4  1  77, 2.  042  57,  G.  ./4  3  3  c,  1 2  N2i 
GGEC  I  STA  =  . 30, 1  ST  e*30,  THK  =  0. 0  1000,  0.  0350C,0.0J5oC,  7  DLXsO.  06 15  1,0.  06059, 0.05  8  96,  G.  O':  7  34,0.0557  1  ,tlM 
f-f'Fo  ISTA  =  31,  ISTB=3  1 ,  THK  =  C.  0 1 000 , 0. 03500, 0. 0  3  500,  7  0LX  =  O.  04055, 0. 0**  100, 0. 04  1  79 , 0.  U4  2  57, »...  v4  j  36 , 1  £NL 
6GEU  ISTA  =  3  2,ISTB  =  32,THK=O.C1000,0.03SOC,0.0  350C,7DLX  =  0.06110,C.  06  026,0.  05b80,O.C5734,C..  0  55  6r;,?  ENL 
LG  SC  IS7A  =  33,IS7B  =  33,'rHK=C.01COO,C.C35yO,O.C3500,7  i:LX=0.040bO,0.04lOU,0.04l  81,0.  04257,O.C4334,G£:ir 
r.GiO  :bTA  =  34,  ISTB  =  3  4  ,THK  =  G.01000,0.03500,G.0  3SOC,  7  OLX  =0.06055,0.  G59c3 ,  G.  v  585s  ,  C ,  057  3  4 , 0. 05t  o9 ,  £  £  NO 
EGEO  ISTA=3S,  1ST  P  =  1S,THK=0.01000,  0. 03500, 0 ,0  3  500,  T  CLX  =  0.  6  ,  G.  C4  1 C9  ,  C  -  04  1  8  3 ,0 , 042  57, 0-  04  2  3  2 , 6  t  .S  : 

r.OF.O  I.STA  =  3  6,ISTo  =  36,THK  =  C.01000,  0.03  500,0. 0J50C,TDLX  =  0.Ce  CGC,  0.  05941,0.  C  58  3  7,0. 057  3i,C.s^6c  3  1,  i  iSO 
EGSD  IS iA=J7, ISTS=37,THK=0. 01000, 0. 03500, 0 . C350C, TCLX=0.0407 3 , C. 04114, C-Ci1ct,C. 04257, c.:43:9,t£N: 
EGEii  :srA=38,  lbT8  =  3H,  'i  hK=G.  0 1 OCO ,  0.  03500, 0. 0  35v7  0,  7  LLX  =  0. 0  5  9^6 ,  «.  0‘-9  ,7  ,  C.  0  ’  8  2C  ,  0 . 0  =  7  ja ,  o.  C  5646,  o  r  N  : 
SGEO  I.STA  =  39,ISTS=3i,IHK  =  C.01000,  0.  0350  0,1).  DJf)00,iDlX=0.C4  CoG,  0.04  119,0.  Cl  Ur,:.  04.5  7,0.  C  4.- 2  7,  £f.S., 
f.G£0  IGTA  =  4  0,  ISTE=4  7,7  HK  =  0.  01000,  0.03500,  0.1  3613,  7CLX  =  0.05j59,0.  05321, 0.05255,0.  05169,  3.  0  5  1  5  3 ,  £  t  N  : 
SGEO  ISTA=4l,  IS7B=4  1 , 1HK=0.0 1  000, 0. 03500,0 . 1 36 13, 7  CLX=0.  0«352, C.  04 392 , 0..  C  44  62 , 0.  C4 53  3 , 0 .  C  4{.  j7 ,  f  E  S L 
EGEO  XSTA  =  4  2, 1ST  8=4  2, 7  HK  =  C.  0 1  000,  0.  03500, 0 . 12  301  ,TDLX=0.  051  1 3 ,0.  05  £82 ,  C.  05027 , 0.  C4  97  3 ,  C .  4  74  4 ,  i.  r  SO 
EGEO  :SIA  =  4  },  157  8*4  3,  THK=0.  01000,  0.03500,0.  12  3C 1 , 1  DLX=0. 04  36  1  ,C.  C4  399,0.  0  44  to ,  0 .  04  f  3  3 ,  U.  :-4  76  4,  t£M, 
EGEO  rSTA=4  4,ISTB=44 ,7HK  =  0.01000,0.03500,‘’.  1 1 038, TCLX=0. C50 1 6 , 0. 04969 , 0. C 4y 42 , G. Cwc 95, 0. C 4 7. c , £  £  N t 
SGEO  1 3T  A*  4  5, 1ST  B  =  4  5,  THK  =  0.0 1  000, 0.03500,0. 1 1 038,TDLX  =  0.  G4370,  C.  04406 , 0-  04470,0.  045j  3,  c-G47^d,tL|,L 
<■0  80  ISTA  =  46, ISTB=46,THK  =  0.0 1000,0. 03  500,0.099 25,7 DLX  =  0. 04962,0. C4 936 , 0- C4 890 , 0 . C4 044  , C . 04 L99 , U  S. 
EGEO  ISTA  =  47, IST£=4 7, THK  =  0.01000,  0.0350C, 0.09825, 7CLX=0. 04379, 0.044  13,0..  04473, 0.04533, f. 04 1 99, fESr 
EGPO  ISTA  =  48,ICTb=4ti,THK=0.01000,  0.03498,0. 08648,7  DLX  =  0..  04923, 0.  04096, 0..  04849,0. 04802, 0.04  6  74,  EENC 
EGEO  ISTA  =  49,I3TB=49,THK=0. 01000, 0.03500,0.08648,1  CLX  =  0. 04388,0. C4420, C«  044 76 , 0. 0453 3, 0..  04674 , EL NP 
SGEO  IS7A  =  5  0,  rSTS=50  ,  T  HK  =0.  0  1000,  0..  03500, 0 . 07494,  1  DLX=0. 049  84 ,0.  04857 , 0..  Cm  b  1 1 , 0.  C47b  4 ,  o.  u4 1 G  1 , 6  i  s: 
SGEO  ISTA  =  5 1, IS73  =  5 1 ,7HK  =  0. 01000,0. 03500,0 . 07a94, 7 DLX  =  0.0 4396 ,C.  04427 , 0 . C4480 ,0. Ou 5 3 3 , 0. 04i.  3  1 , 1  I v T 
SGEO  ISTA*  52, 1ST  B  =  62,THK  =  C.  01  000, 0.03  500,0 .0636  1,7  DLX  =  0.  04836,  0.  048  1 0,  0 .  G4  764  ,  C .  04?  1  e  ,  C.  G  4t  2y  ,  i  E  .i: 
LGEO  ISTA=53,IST6*53,THK=O.0lCOO,O.O35OO,0.O63bl,7tLX=O.O44O5,O.C4433,O. OU463,0.0U5i3,C.04c29,f EM 
SGEO  ISTA  =  54,ISTB=54,THK-0.010C1,O.C  3  502,  ).  C52  38, 1  DLX  =  0.  04794  ,  C.  04  7b9  ,  C..  C4  724 , 0 .  04 1  c  C,  S  .  G4  r  ^  o ,  I  E  {.  r 
SGEO  ISTA=55, 1ST E=55,THK= 0.0 1000, 0.03499,u.0523e,lLLX=0.04412,0. 04U38,Q.044e6,C.C4533,j.04bOb,S£M 
SGEO  ISTA=56,ISTB=5b,THK=0. 01007, 0- 035 1 0,0 . 04 108, TCLX=0. 04868 ,C. C4e39 , G. 04 7 c4 , 0. 0m7 It , 0. j 46 1 ? , f £ sT 
SGEO  ISTA=57, ISTE=57,THK=0. 00999,0.03496,0. 04 108,TDLX=0. 04400,0. 04426,0.04477,0. C453 3, 0. . 4 u 1 7 , E i NU 
SGEO  ISTA=58,ISTB=5  6,THK=0.00915,  0.  C3596 ,0 . 03001 ,  TDLX  =  0. 0446  1 , 0,  04468 , 0.  04  5  1 9  ,  G.  C  u '■yb,  i  .  3  46  .  i ,  t  E  St 
SGEO  ISTA  =  59,ISTB=59,THK  =  0-01001, 0.0350 3 ,0 . 0300 1 , TCLX  =  0. 04250 ,  C.  04 2 8 3 , 0- 04  j c2  ,  c  .  045 i  j,  w.  c  4L2 *• ,  t  h  Si: 
SGEO  ISTA=60, ISTB=60, T HK  =  0.00728,  0.03426, 0.03C04,7CLX  =  0,04 119,0.  04100,0.0406 7, j.C4 11. ',..041  OS, t EM 
SGEO  ISTA*61,ISaB=61 ,THK  =  0.008 13, 0-030 77, 3.0 "004,7 CLX=0. 037 8 8,0,  03819,0.039  13,0.  0403 1 , 0. 04  1 05, E £  ST 
SGEO  ISTA  =  6  2,ISTB=62, THK  =  0,00  551, C-031  72 , 0. 02  998, 7DLX=0, 04080, 0.  04Cb9,  0.  04  045,0. 04036,0. 040 3 7, EE  NO 
SGSO  ISTA=63,ISTB=63,THK=0. 006 18,0. 02736, 0.C2998,1DLX=0. 03978, 0.03993, 0.04 01 7, 0.04 03 1,0.U4wJ7,fctM; 
SGEO  ISTA =64, ISTB=64,THK=0- 00375,0.02846,0. 02998,1 DLX=0,040o 5, 0. 04057,0.04037,0. 04C31 , 0.G4O3 1 , FEND 
SGEO  ISTA  =  65, 1ST  B=65,THK  =  0.  004  19, 0-02  459,0.  C2598,TDLX  =  0.- 04008, 0.0401 5, 0.04028,0.  0403  1,0.  04  03  1  ,  EEM 
SGEO  ISTA=66,ISTB=66,THK=0- 00  198, 0-02452, 0.  (2998, 1 tLX=0. 04058 ,0-  04 053 , 0. 04037 , 0. 0403 1 , 0. 040 3 1 , 6 1 N.. 
SGEO  I STA=67, 1STB=6 7,THK=0.002 18, 0.022 19, 0. 02 S98, 7 CLX  =  0.0402 4,0. 04 C 25,0-0  4030,0.  04031,0-040 11, ft  NT 
SGEO  ISTA  =  68,ISTB=6  8,THK=O-OOO23,O.O1991,O.O2S90,TDLX=O,O4O4  9,O.O4C4b,O.O4O35,O.O4O31,:.  J4:3  1,£lNC 
SGEO  I STA  =  69, 1 STB=69,THK=0. 000 18, 0.01978,0. 02998,1 CLX=0. 0403 5,0. C4C35, G.  04  032,0.  C4C3 1 , j .  04C  M  ,  Sr  ;o 

(f)  TACTl  input  for  blade  slice. 

Figure  13.  -  Continued. 
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Figure  13.  -  Concluded, 
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Figure  21.  -  Calling  relation  of  subroutines.  Called  subroutines  are  always  below  the  calling  subroutine.  (This  is  not  a  flow  chart 
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